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M.  ABSTRACT  ^ComMwo  on  r'«^oo  nt4m  II  nocooMry  mmthmt) 

Physical  characteristics  of  the  head  and  neck  were  measured  on  18  male 
Navy  volunteers  who  had  previously  undergone  testing  on  the  NAJ^RL  sled 
facility  at  Michoud  Station,  New  Orleans.  Measurements  include  55  standard 
anthropometric  measures,  32  anthropomet. ic  measures  of  the  seated  subject, 
three  dimensional  head  and  neck  range  of  motion,  neck  muscle  reflex 
times  in  response  to  head  jerks,  and  neck  muscle  voluntary  isometric 
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strength.  These  latter  measurenents  were  tai:en  in  both  the  sagittal  and 
lateral  planes. 

heasure'jent  results  were  used  to  establish  paraaeter  values  for  the 
NVMA'2D  Crash  Victim  Simulator  data  set  in  an  attempt  to  reproduce  the 
dynamic  response  of  these  volunteers  to  sled  accaleratio.i  at  6 and 
15  G's.  Procedures  used  for  computing  the  various  puraneter  values  ord 
comparisons  between  predicted  and  experimental  results  are  presented. 

In  addition,  measurement  data  for  18-24  year  females  taken  p/^iously  have 
been  utilized  to  predict  the  dynamic  response  that  would  be  expected  if 
tliese  subjects  were  tested  at  6 and  15  G's. 

rurther  work  in  studying  the  significance  and  relation  of  various 
physiological  and  biomechanical  parameters  and  of  stimulus  and  experi- 
mental test  conditions  to  the  dynamic  response  is  planned  using  both 
modeling  and  correlation  techniques.  Measurement  data  for  oth^  segments 
of  the  adult  population  will  be  used  to  extetid  the  NANRi.  results  to  the 
general  adult  occupant  population. 
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SUMMARY 


Physical  characteristics  of  the  head  and  neck  were 
measured  on  18  yoriig  adult  male  Navy  volunteers  who  had 
previously  undergone  tests  on  the  NAMRL  sled  facility 
in  New  Orleans.  Measurements  taken  include  55  standard 
anthropometric  measures,  32  cmthropomct.;ic  measures  of 
the  seated  subject,  three  dimensional  head  and  neck  remgc 
of  motion,  neck  muscles  reflex  tiroes  in  response  to  head 
jerks,  and  neck  muscle  voluntary  isometric  strength.  The 
reflex  time  and  strength  tests  were  performed  in  both  the 
sagittal  and  lateral  planes.  The  ramge  of  motion  results 
for  this  gro\ip  of  18  NAMRL  subjects  were  in  good  agree- 
ment with  results  for  18-24  year  males  and  females  from 
the  general  population.  In  the  sagittal  plane,  the 
a/erage  range  of  motion  angles  in  extension  were  79.0  and 
60.5  degrees  respectively  as  measured  from  the  Frankfort 
Pl2me  position.  Reflex  times  were  similar  for  flexion, 
extension,  and  lateral  bend,  being  53.5,  55.5,  and  51.5 
msec  respectively.  In  strength,  the  group  of  NAMRL 
s’ibjectu  was  similar  to  35-44  year  males  of  the  general 
population.  The  greatest  strengths  were  in  extension 
where  the  average  is  about  33%  greater  than  in  flexion 
or  Jateral  bend. 

Where  appropriate,  these  measurement  results  for  5 of 
the  18  subjects  were  utilized  in  establishing  a data  set 
for  the  MVMA-2D  Crash  Victim  simulator.  Simulations  of 
NAMRL  sled  te.^ts  at  -Gx  impact  accelerations  of  6 and  15 
G's  were  made  using  either  the  experimental  T^  accelera- 
tions as  input  to  the  neck  or  experimental  sled  accelera- 
tion profiles  as  input  to  the  sled.  Simulation  results 
for  head  angula'  acceleration,  head  angular  velocity,  head 
angular  position,  head  resultant  acceleration,  and  Tj 
resultant  acceleration  are  compared  with  average  experi- 
mental results  out  to  300  msec  for  the  group  of  five 
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SiU:.jects.  Results  to  date  indicate  reasonably  goou 
agreeinent  between  experimental  and  simulation  curves 
at  both  6 and  15  G's-  Further  work  is  needed,  however, 
to  improve  certain  aspects  of  the  model  such  as  joint 
stop  characteristics,  passive  tissue  modeling,  2md 
restraint  system  iDodeling.  Effects  of  changing  muscle 
tension,  chest  compli2mce,  joint  stop  stiffness  co- 
efficients, and  upper  torso  joint  stiffness  (i.e.,  amount 
of  torso  flexion)  have  also  been  examined  using  the  MVMA- 
2D  model.  Results  obtained  with  varying  aunounts  of 
muscle  tension  indicate  that  muscle  effects  are  more  pre- 
dominant at  6 G's  thcin  at  15  G's. 

Simulations  for  18-24  year  females  at  6 and  15  G's 
were  made  using  measurement  data  obtained  in  previous 
studies  at  KSRI.  Results  were  not  dramatically  different 
from  the  NAMRL  simulations,  the  primary  difference  being 
an  increase  in  the  maximum  flexion  angle  of  the  head  by 
about  40  and  25  _ ercent  at  6 and  15  G's  respectively. 


Chapter  1 

INTRODUCTION  AND  OBJECTIVES 

A.  Statement  of  Project  Goal 

Measurements  of  dynamic  responses  to  impact  acceler- 
ations which  have  been  taken  on  a selected  male  military 
population  at  the  Naval  Aerospace  Medical  Research  Labora- 
tory (NAHRL)  at  Michoud  Station,  New  Orleans,  represent 
the  most  cotaprehensive  source  of  information  availeU:;le 
related  to  the  dynamic  response  of  the  human  head  and 
neck.  To  what  extent  these  data  represent  the  total  adult 
U.S.  population  is  unknown,  however.  In  recent  studies 
sponsored  by  the  Insurance  Institute  for  Highi’ay  Safety 
(IIHS)  and  conducted  by  the  Highway  Safety  Research  Insti- 
tute (HSRI) , basic  information  which  is  believed  to  be  re- 
presentative of  neck  physical  Ghjiracteristica  for  the  adult 
U.S.  jxjpulation  from  18  to  75  years  has  been  obtained. 
Included  in  thesu  data  are  anthropometry,  head/neck  range 
of  motion,  neck  muscle  strength,  and  neck  muscle  reflex 
time  measurements.  The  primary  purpose  of  this  study  is 
to  determine  to  what  extent  these  data  may  be  used  with 
mathematical  modeling  techniques  in  order  to  extend  and 
project  the  NAMRL  dynamic  respon.se  results  io  the  general 
adult  U.S.  population. 

B . Background 

Response  of  the  human  head  and  neck  to  imf>act 
accelerations  is  a matter  of  major  concern  in  the  design 

The  riglits,  welfare  and  informed  consent 'oi?  t^ie 
volunteer  subjects  who  participated  in  this  study  were 
observed  under  guidelines  estedslished  by  the  U.S. 

Department  of  Health,  Education  and  Welfare  policy  on 
protection  of  human  subjects  and  accomplished  under 
medical  research  design  protocol  standards  approved  by 
the  Committee  to  Review  Grants  for  Clinical  Research  and 
Investigation  Involving  Human  Beings,  Medical  School, 

The  University  of  Michigan. 
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of  biomechanical  models,  anthropomorphic  duimies,  and 
occupant  crash  protection  devices. 

There  are  a large  number  of  studies  which  have 
attempted  to  determine  the  relationships  of  head  injury 
and  concussion  to  impact  forces,  but  only  recently  has 
attention  been  given  to  the  respective  influences  of  the 
effects  of  head  motion  upon  injury.  It  is  still  unclear 
what  the  respective  erfects  of  rotational  and  translational 
forces  may  be.  Results  of  experiments  reported  by 
researchers  such  as  Holbourn  (14),  Pudenz,  et  al.  (19), 
Martinez  (16) , and  Ommaya  (18)  have  indicated  that  rotation 
alone  can  cause  brain  injury  and  concussion  in  whiplash. 
However,  Hodgson  (13)  , Gurdjian  (12)  and  others  contend 
tiiat  other  factors  such  as  resultant  intracranial  pressure 
gradients  may  cause  trauma  by  high  shear  stress  concentra- 
tion in  the  brain  stem  and  upper  spinal  cord.  Young,  et  al. 
(26)  have  recently  demonstrated  concussion  to  the  fixed 
primate  head  without  translational  movement.  Clarke  (4) , 
studying  human  volunteers  in  dynamic  tests  of  adult  males 
at  peak  sled  velocities  of  26.2  ft/sec  and  V.8  to  10  G, 
concluded  that  peak  head  angular  accelerations  and  linear 
resultants  may  have  less  trauioatic  consequences  than  the 
degree  of  head-neck  hyperex tens ion. 

This  disagreement  among  researchers  as  to  the  mecha- 
nisms of  injury  in  head  impact  and  whiplash  is  also  seen 
in  consideration  of  critical  values  of  rotational  velocity 
and  acceleration  at  which  concussion  occurs  in  man.  Recent 
work  by  Ewing  and  Thomas  (6)  using  male  human  volunteers 
in  dynamic  sled  tests  found  no  clinically  observable  effects, 
due  to  acceleration  on  a subject  in  which  the  peak  mouth 
angular  velocity  exceeded  30  rad/sec  (at  10  G,  250  G/sec) , 
although  this  level  had  been  previously  considered  by 
M.hone,  et  al.  (15)  and  Ommaya  to  be  the  critical  level 
for  human  concussion. 
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similarly,  although  there  has  been  considerable 
effort  to  realistically  simulate  the  human  neck  in 
various  versions  of  an  "improved"  anthropomorpnic  dunmiy, 
the  lack  of  valid  human  bioengineering  data  has  remained 
a major  problem  and  much  controversy  in  this  area 
continues. 

Thus  the  continuing  series  of  . npact  acceleration 
tcKts  being  conducted  by  Ewing,  et  al.  (5-10)  using  human 
voliiitoor  subjects  have  been  of  particular  significance 
since  this  effort  has  resulted  in  an  extensive  body  of 
kinematic  experimental  data  under  dynamic  conditions. 

This  wi..rk,  which  has  involved  precise  measurement  of  the 
complete  input  accelaration  to  the  head  and  neck  (measured 
at  the  first  thoracic  vertebra)  , precise  naasurament  ‘';f 
the  dynamic  response  of  the  head  and  neck  tc  he  inpv\t 
accel»ration , and  development  of  dat.c?  acqoisicion  and 
automatic  processing  fl^sveiris,  must  be  characteri;:»vd  a« 
producing  the  most  extensive  dynamic  data  vising  the  most 
sophisticated  experimental  techniques  and  precise  instru- 
mentation to  date  for  the  impact  range  under  stud/. 

Primary  objectives  of  the  NAMRL  research  effort  arc 
to  acquire  data  that  can  be  used  to  1)  develop  design 
criteria  for  construction  of  dvimmies  which  will  closely 
reproduce  man's  response  to  crash  acceleration  , anr’ 

2)  define  the  envelopes  of  impact  acceleration  which 
result  in  the  injury.  If  these  results  are  to  include 
concern  for  the  total  population  who  may  be  involved  in 
crash  situations  in  both  military  and  non-milit«a;y  vehicles, 
then  it  becomes  important  to  be  able  to  extend  these 
dyneunic  response  data  to  the  general  U.S.  adult  population. 

In  two  studies  by  Snyder,  et  al.  (22,23)  and  rej^ortod 
by  Foust,  et  al.  (11)  and  Schneider,  et  al . ^■^l)  , basic 

data  concerning  physical  characteristics  of  the  head  and 
neck  were  obtained  on  a sample  of  subjects  designed  to 
represent  the  adult  vehicle  occupant  population.  Those 
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studies  were,  to  our  knowledge « the  most  coaprehensive 
attempt  to  relate  such  physical  measurements  such  as 
muscle  isometric  strength , muscle  reflex  time,  cervical 
range  of  motion,  and  cuithropometry  to  the  age,  sex,  and 
stature  of  a population  representative  of  U.S.  adults. 

While  subjects  vere  tested  in  both  sagittal  and 
lateral  planes,  the  general  relationships  of  the  measured 
physical  characteristics  to  age  euid  sex  were  the  same  in 
both.  Cervical  range  of  motion  was  greatest  in  the 
rotational  plane  and  smallest  in  the  lateral  plane,  and 
shovred  an  average  decrease  with  age  of  20-45  percent  from 
young  to  elderly  subjects.  Neck  muscle  reflex  times 
rangevi  from  eibout  30  to  75  msec,  were  generally  smaller  for 
lateral  head  movements,  showed  increase  with  sxibject  age, 
and  were  slightly  shorter  on  the  average  for  fescales. 

Muscle  strength  was  found  to  be  about  33%  greater  in  ex- 
tension than  in  flexion  or  lateral  bend,  showed  a decrease 
with  age,  and  was  on  the  f/erage  1-1/2  to  2 times  greater 
in  males  than  in  fonales.  No  significant  correlations 
oet\/fejr.  these  measurements  and  subject  anthropometric 
lauus’tres  were  found. 

A basic  assumption  underlying  the  application  of 
these  data  to  studias  o*j  hu'u^ul  impact  tolerance  is  that 
there  is  a relationship  between  the  differences  j.i 
physical  characteristics  of  individuals  and  dif j erences 
in  dynauKiic  response  to  impact.  If  this  assumption  is 
valid,  a potentially  productive  r(‘.'^^;ecrch  program  would 
involve  bridging  the  gap  between  the  dynamic  studies  of 
a highly  selected  population  on  the  one  hand,  and  the 
essentially  static  rceasurecaents  representative  of  the 
U.S.  adult  population  on  the  other.  The  present  study 
was  undertaken  out  of  these  considerations. 
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C.  Objectives 
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In  order  to  accomplish  the  goals  of  this  study, 
two  principal  objectives  must  be  achieved.  First,  it 
must  be  shown  that  the  static  measurements  can  be  utilized 
in  a mathematical  model  to  give  accurate  simulations  of 
experimental  impact  results,  and  that  this  validated  model 
can  be  used  across  a range  of  impact  accelerations  and 
conditions.  Secondly,  it  must  be  demonstrated  that 
reason'jijle  correlations  exist  between  measured  physical 
properties  and  experimental  dyn^unlc  response  characteristics. 
In  order  to  attain  these  objectives.-  the  following  tasks 
were  established. 

1)  Conduct  all  of  the  HSRI  sagittal  plane  and 
lateral  plane  testing  on  a group  of  NAMRL  subjects  who 
have  previously  undergone  acceleration  impact  testing 
over  a range  of  acceleration  levels.  iQata  collected 
would  include  standard  anthropometry,  seated  anthropo- 
metry, three  dimensional  voluntary  range  of  motion  of 
the  head  and  neck,  neck  muscle  stretch  reflex  times  and 
acceleration  in  response  to  head  jerk,  and  n'ayimum 
voluntary  neck  muscle  isometric  strength. 

2)  Compare  NAMRL  measurement  results  •;  reJuTts 
from  the  O.S.  adult  population  obtained  in  the  1. 
sagittal  and  lateral  plane  studies  (11,  21). 

5)  Use  the  NAMRL  measxireraent  results  to  establish 
parameter  values  where  appropriate  in  the  data  set  of 
the  MVMA-2D  Crash  Victim  Simulator. 

4)  Run  the  MVMA-2D  model  with  this  data  set, 
appropriate  initial  conditions,  and  stimulus  inputs  in 
an  attempt  to  reproduce  the  experimental  results  from 
NAMRL  subjects  at  several  acceleration  levels.  Adjust 
or  "tune"  other  parauneters  for  which  no  data  are  available 
in  order  to  obtain  optimal  matching  to  experimental  curves. 
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5.  Determine  if  and  %^en  any  relationships  exist 
between  static  measurements  and  dynamic  response,  and 
confirm  that  these  relationships  are  supported  by  the 
model  by: 

a)  correlating  various  static  measurements 

with  various  peak  pareuaeter  valuej  in  experimental 
response  curves  for  the  group  of  SAMRL  subjects 
measured. 

b)  eoLamining  changes  in  response  curves  for 
subjects  with  different  physical  characteristics 
and  using  these  different  measurements  in  the 
model  to  see  if  the  changes  are  predicted. 

6.  Use  the  measurements  obtained  in  the  IlHS 
studies  to  predict  the  dynamic  response  results  that 
would  be  obtained  if  other  segments  of  the  population 
were  tested  at  the  NAKffL  sled  facility. 

7.  Use  the  validated  MVMA-2D  model  to  predict  the 
response  of  occupants  subjected  to  more  realistic  axxi 
practical  crash  situations. 

8.  Use  the  validated  KVMA-2D  sK)del  to  predict 
the  response  of  NAMKL  subjects  to  sled  tests  where  the 
acceleration  vector  is  of  a greater  magnitude  than  can 
be  safely  used  with  volunteer  subjects. 

At  the  time  of  this  report,  tasks  1 through  3 have 
been  completed  for  a group  of  18  NAMRL  subjects  who 
have  undergone  sled  tests  up  to  15  G's  in  the  -Gx 
direction.  Considerable  progress  has  been  made  with 
task  4 although  further  work  and  improvements  in  the 
model  are  needed.  Measurement  procedures  and  results 
are  presented  in  Chapter  2 while  Chapter  3 describes  the 
procedures  emd  results  to  date  concerned  with  tasks 
3 and  4.  Chapter  4 gives  the  results  obtained  by  using 
the  IIHS  data  to  predict  sled  test  responses  at  6 and 


15  G'e  for  18-24  year  fenjiles  (tesk  6).  Chapter  5 containa 
a brief  discussion  of  the  results  to  date  and  suggestions 
for  future  work  on  tasks  4 through  8. 


Chapter  2 


MEASUREMEHT  OP  P , :CXCAL  CHARACTERISTICS 
A.  General 

Liglitoen  male  N:i  / personnel  who  had  previously 
undergone  testing  cn  the  NAMRL  sled  facility  were  brought 
to  the  Highway  Safety  Resesirch  Institute  (HSRI)  for  two  to 
three  days  for  measurements  of  physical  characteristics 
related  to  the  head  and  neck.  The  subjects  were  brought 
to  HSRI  in  three  groups  of  7,  6,  and  5 over  a period  of 
one  month.  They  were  assigned  subject  identification 
numbers  consisting  of  a prefix  code  denoting  the  sex,  age, 
and  stature  of  the  individual  and  a chronological  testing 
number  from  1 to  18.  Table  2.1  is  a list  of  these  HSRI 
subject  codes  together  with  the  corresponding  NAMRL  subject 
number.  The  letter  N denotes  NAMRL  male  subject  (to  dis- 
tinguish from  others  previously  tested) ; A indicates  that 
the  subject  was  between  18-24  years;  while  S,  M,  and  T 
correspond  to  short,  medium,  and  tall  according  to  the 
l-20th,  40-60th,  and  80-100th  percentiles  of  stature  for 
this  age  and  sex. 

While  the  immediate  concern  in  this  study  was  to 
measure  those  physical  characteristics  which  relate  to 
head  2Uid  neck  movement  in  the  sagittal  plane  (since  NAMRL 
sled  testing  is  with  -Gx  acceleration) , the  capability 
existed  from  pre\lcus  studies  to  test  in  lateral  bending 
and  these  measurements  were  also  taken  on  each  subject  for 
future  use.  Measurements  taken  include  standard  anthro- 
pometry, anthropometry  of  the  seated  subject,  head  and 
neck  range  of  motion,  neck  muscle  stretch  reflex  times, 

2md  neck  muscle  isometric  strength  capability.  The 
following  s€K:tions  describe  briefly  the  testing  procedures 
used  (these  are  discussed  in  detail  in  references  22 
and  23)  and  also  present  the  measurement  results  for  this 
group  of  18  Navy  subjects. 
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TABLE  2.1 


HSRI  A.:;.  luAMRL  SUBJECT  NUMBERS  FOR  18  NAVY  SUBJECTS 


HSRI 

.\AMRL 

NAMCl 

K-39 

HAT03 

]i-l9 

H-h? 

!.'A,VG5 

NAVOr'' 

:;asc7 

H-33 

H-32 

NAMG9 

4. 

:;AMir 

H— 7 

■ ml  ' 
«iM  A X ▲ 

H-^1 

:.A'C.r 

H-50 

NAMl  ^ 

HAT  I-'. 

H— 1 

:,a:’:  ^ 

M _ *?  J 
' • J / 

r.'AMit' 

K-33 

nat:t 

H— ^0 

\AT16 

H-3T 

Pramliq;  paft  Haik 
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B.  Anthropometry 


1.  Methods . A total  of  87  anthropometric  measure- 
ments were  obtained  on  each  subject  during  the  initial 
phase  of  testing.  Figure  2.1  is  a list  of  these  measure- 
ments divided  into  two  groups.  Group  I contains  55 
measurements  taken  by  standard  techniques  to  describe  the 
general  bod^  characteristi  the  head  and  neck,  location 
and  sizes  of  body  masses,  and  body  somatotypes.  Group  II 
contains  32  measurements  taken  to  describe  the  position  of 
the  seated  occupant.  Of  these  32,  7 were  taken  with 
standard  equipment  while  25  are  measures  tcLken  or  derived 
from  orthogonal  photogrammetry  techniques.  The  photo- 
grammetry  setup  consists  of  a set  of  th,  ee  Pentax  cameras 
orienced  orthogonally  to  each  other  and  aimed  toward  the 
subject  from  the  front,  left  side,  and  top.  The  cameras 
are  aligned  such  that  the  centers  of  their  focal  planes 
intersect  at  a common  origin  within  the  subject.  The 
locations  of  high  contrast  markers  placed  on  the  subject's 
head  and  upper  torso  and  visible  in  at  least  two  cameras, 
may  be  determined  in  three  dimensions  by  geometric 
relations  and  measurements  from  the  films.  The  films  are 
projected  onto  a tablet  digitizer  and  the  points  are 
digitized  in  a specified  sequence  onto  paper  tape.  The 
tapes  are  later  analyzed  by  computer  programs  which  com- 
pute the  3-dimensional  location  of  the  points. 

2.  Results . Table  2.2  presents  the  statistics  of 
the  anthropometric  measurements  for  the  18  NAWRL  subjects 
measured  while  Table  2.3  gives  the  statistics  for  the  three- 
space  location  of  the  head  and  torso  of  the  seated  subject. 
In  many  cases,  the  names  of  measurements  have  been  abbre- 
viated for  compactness  and  these  abbreviations  may  ha 
cross  referenced  with  a more  complete  measurement  name  in 
Appendix  A.  Tables  B.l  and  B.2  in  Appendix  B contain  the 
individual  data  for  each  subject  from  vhich  these  statis- 
tics were  derived. 
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TABLE  2.2  ANTHROPOJBETRY  STATISTICS 
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kl  IKC) 
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UB) 
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SIAl  (C£i 

Id 

pcaiihix 

Id 

EaSllt'i 

18 

tlkLClB 

Id 

tikLlLib 

Id 

EITEGDI 

Id 

Id 

bkkLLQ 

Id 

BA'^AbL 

17 

CCetAhC 

Id 

alluGGLE 

1b 

fliificetii 

1b 

bllhQllk 

18 

f ACi  fcl 

lb 

lAXlKEf 

1b 

IPUKtb 

Id 

SUIlKCIt 

16 

iiiJi.KCIE 

1b 

tCSlkKIC 

1b 

El AC  a c P 

lb 

EIDiIXfca 

16 

CEESJET 

16 

CEESlEii 

18 

CLESICIE 

16 

WAISiEI 

1b 

iiAISltb 

16 

•‘A1B1QU 

1b 

hlPhi 

1b 
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1b 

H i PCI « 

1b 

ACEhACLG 

1b 

AeHCltAH 

1b 
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ElCf iCiE 

1b 
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1b 
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1b 
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16 
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1.9 

65.6 

2.  4 

13.  S 
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15.0 

C.  4 

19.6 

C.6 

35.  1 
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w • £ 
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1.3 

9c.  6 
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9 fa  . •« 
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132.  2 

3,6 
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2,  1 

95.0 

e.  1 

106.7 

2.8 

29.  S 

3,  3 

b3. 9 

6.3 

92. 7 

2.  7 

33.  t 

2.  C 

96.5 

6.  3 

32.  6 

1.6 

32.9 

3.  3 

2t.  1 

3.  7 

33.  1 

-»  2 

26.  1 

t.  9 

26.  4 

1.6 

17.4 

0 . b 

18.  9 

0.  4 

41.2 

1.  5 

57.5 

5.  S 

39.  3 

3.  7 

40.  6 

1.6 

45.  7 

1.  6 

MiNlEU.l 

MAXIi'JOfl 

61.1 

105.  2 

134.5 

2 31.  5 

165.5 

1d4.  4 

23,7 

34,  9 

65.0 

9 9.8 

54.6 

b 1.  3 

62.2 

70.  5 

U.  J 

14.  d 

14.4 

15.  9 

13.0 

20.8 

32.  1 

37.  6 

33.0 

Jo.  7 

27.7 

3 1.  3 
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34.  5 

24.8 

32.  3 
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1 4 . 0 

11.3 

13.  3 

10,2 

I4.  9 

34.9 

44.  3 

35.6 

4 3.  9 

l4. 9 

7 9.  5 
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44,  9 

4J,3 

54.  2 
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Ud.  3 
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34.  7 
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26.5 
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4 1 , 6 

4 9,6 

16 


TABLE  2.2  (continued) 
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T'PiuE  2.3  UPPm  TORSO  AND  HEAD  LANDHASK 
COCOiDlSATE  STATISTICS  (re  SRP) 
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C . Range  of  Motion 


1.  Methods.  As  described  in  references  21  and  22, 
orthogonal  photogranscetry  was  also  used  to  determine  the 
subject's  head  and  neck  range  of  motion.  Immediately  fol- 
lowing the  seated  anthropometry  measurements  the  subjects 
were  asked  to  perform  the  sequence  of  head  ntovenents  shown 
in  Table  2.4.  In  each  position  when  the  subject  attained 
the  limit  of  his  voluntary  movement,  photographs  were 
taken  simultaneously  by  the  three  cameras.  Figure  2.2 
shows  a subject  performing  the  head  and  neck  extension 
movement.  The  Eulor  angles  describing  each  position 


Figure  2.2  NAMRL  Subject  Performing  Range-of-Motion  Tests 


relative  to  the  Frankfort  plane  position  were  computes^  by 
digitizing  the  points  on  the  coordinate  syatean  head  piece 
worn  by  f e subject  throughout  these  movements.  The  vec- 
tors desc7 ibing  the  orientat  on  of  these  coordinate  axes  in 
opacfc  were  determined  in  a manner  similar  to  the  seated  an- 
thrcj’xxscatric  measures  and  the  Euler  angles  computed  by  ap- 
propr.'.ate  ev'uations  (see  appendix  E of  reference  22  ) . 
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TABLE  2.4 


SEQUENCE  OF  RANGE-OF-HOTION  POSITIONS 

1 . Frankfort  Plhuie 

2.  Normal 

3.  Extension 
U.  Flexion 

5.  Right  Rotation 

6.  Left  Rotation 

7.  Right  Latersd  Bend 
6.  Left  Lateral  Bend 

9.  Left  Rotation  Plus  Bend  Toward  Lert 

Id.  Left  Rotation  Plus  Bend  Tovard  Rear 

11.  Right  RotatJor  Plus  Bend  Toward  Left 
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In  computing  the  Euler  angles,  the  order  of  movements  to 
attain  a given  position  is  asnunved  to  be  yaw  (rotation)  , 
pitch  (flexion  or  extension) , and  then  roll  (lateral 
bend)  and  the  Euler  angle  reference  frame  2oces  are  as 
shown  in  Figure  2.3  where  positive  x is  forward,  positive 
y is  toward  the  right,  and  positive  z is  down. 


K 


K 


Figure  2.3  Euler  Angle  Reference  Frame  and  Angle 
Directions 


2.  Results.  While  the  measures  of  sagittal  plane 
range  of  motion  are  of  primary  significance  to  this  study, 
presentation  of  all  the  results  is  given  here  for  com- 
pleteness and  future  reference.  Table  2.5  gives  the 
statistical  summary  for  these  18  subjects  showing  the  mean 
standard  deviation,  and  minimum  and  maximum  value  for 
each  Euler  angle  at  each  position.  The  significant  num- 
ber is  the  value  for  the  ang''e  related  to  the  plane  of 
primary  movement  (e.g.,  the  pitch  angle  in  flexion)  but 
the  other  angles  give  an  indication  of  the  deviations 
from  this  pl^me  which  might  be  caused  by  performing  the 
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RANGE-CaF'-MOTION  Sl'ATISTICS 
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movement  incorrectly  or  forcing  against  the  physiological 
"stops."  As  with  the  anthropometry,  the  abbreviated  neunes 
can  be  cross  referenced  with  the  list  in  Appendix  A.  Table 
2.6  is  an  attempt  to  simplify  these  results  and  shows  only 
the  average  Euler  angles  for  the  group.  Individual  range- 
of-motion  results  cam  be  found  In  Table  B.3  of  Appendix  B. 

D.  Reflex  Times  and  Strength 

1 . Methods. 

a.  Reflex  time.  Neck  muscle  reflex  times  were 
measured  by  recording  both  head  acceleration  and  neck 
muscle  electromyograph  (EMG)  signals  in  response  to  a 
head  jerk  produced  by  dropping  a 1 lb.  weight  approxi- 
mately 6 to  8 inches.  The  set-up  for  testing  of  the 
extensor  muscles  (splenius  capitus)  is  shown  in  Fitjure 
2.4.  A line  attached  to  a band  placed  about  the  subject's 


Figure  2.4  NAMRL  Subject  Ready  for  Sagittal  Reflex  Time  Test 


head  is  draped,  over  a pulley  and  threaded  through  the  drop 
weight  which  is  held  xn  position  by  an  electromagnet.  When 
a switch  on  a control  console  is  deprefsed,  the  weight  is 
released  and  caught  by  the  small  pretension  or  stop  weight 


TABLE  2.6 


average 

EULER  ANGLES  FOR  18 

SAMRL  SUBJECTS 

AVERAGE  EULER 

ANGLE  RE  FRANKFORT 

POSITION 

PGSTTIOM 

YAW 

PITCH 

ROLL 

Normal 

-.7 

-2.0 

.1 

Extension 

- 

79.0 

- 

Flexion 

1.2 

-60.5 

-2.1 

R.  Rotatlcn 

76.1 

-2.1 

5.9 

L.  Rotation 

-77.5 

-2.7 

-9.5 

R.  LaterxL  Bend 

4.1 

.4 

45.0 

L.  Lateral  Bend 

-4.6 

-.7 

-47.1 

L.  Rot.  + FT.exlon 

-71.3 

-31.1 

-12.7 

L.  Rot.  ♦ LLB 

-73.8 

3.5 

-40.5 

R.  Rot  * Extenaion 

77.0 

23.3 

14.8 
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producing  a “tug"  to  the  sxibject's  head.  Two  surface 
electrodes  placed  over  the  muscle  group  of  interest  (the 
splenius  capitus  for  extensors/  the  stern2imastoids  for 
flexion  euad  lateral  bend)  measure  the  electrical  activity 
resulting  from  the  stretch  reflex  response  of  the  muscle 
produced  by  the  jer)c.  Head  acceleration  was  measured  by  a 
set  of  4 linear  accelerometers  which  were  oriented  in  the 
plane  of  the  head  jerk  in  a configuration  which  allows  for 
calculation  of  resultant  head  angular  and  linear  accelera- 
tions. These  accelerometers  are  mounted  to  a bar  and 
fixed  to  a bite  plate  held  in  the  subject's  mouth  during 
testing.  The  plate  is  fitted  to  each  subject  using  a 
thermoplastic  moldable  dental  compound. 

Each  subject  was  tested  in  sagittal  flexion  (extensor 
muscles) , sagittal  extension  (flexor  muscles) , and  lateral 
flexion  to  the  left.  A series  of  six  or  more  drops  were 
performed  in  each  position  and  the  average  reflex  time 
computed.  Prior  to  each  test  the  subject  was  instructed 
to  relax  and  close  his  eyes,  but  to  attempt  to  maintain 
his  head  erect  when  the  tug  was  felt. 

Figure  2.5  shows  a typical  result  produced  by  the 
weight  drop  where  only  one  accelerometer  signal  is  needed 
and  used  to  compute  the  reflex  time.  The  beginning  of 
muscle  electrical  activity  is  indicated  by  a sharp  spike 
in  the  relaxed  EMG  signal,  followed  by  intermittent 
electrical  activity.  The  time  from  onset  of  head  acceler- 
ation to  this  first  spike  is  called  the  muscle  reflex 
time.  It  is  not,  however,  the  time  required  to  develop 
maximum  muscle  force  which  must  include  a contraction  time 
of  aj^roximately  100  msec.  The  total  time  from  initial 
head  movement  to  maximum  muscle  force  is  therefore  the 
sum  of  the  reflex  time  and  the  contraction  time  and 
could  be  called  a reaction  time. 
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Figure  2.5  Typical  EM6  and  Acceleration  Signals  in 
Response  to  Reflex  Test  Weight  Drop. 


b.  strength.  M&xiiaum  voluntary  isonetric  neck 
muscle  strength  was  measured  on  each  subject  as  a SMuisure 
of  the  strength  capability  of  the  neck  muscles  for  res- 
training  the  head  during  impact.  Tests  were  performed  in 
extension,  flexion,  and  left  and  right  lateral  bend  exer- 
tions with  the  subject  seated  in  the  same  test  seat  as 
used  for  range  of  motion  and  reflex  testing.  Figure  2.6 
shows  a subject  being  tested  for  flexor,  extensor,  and 
lateral  muscle  strengths.  A band  placed  about  the  head  is 
attached  by  an  adjustable  ler«gth  inelastic  rope  to  the 
rigid  test  frenne  via  a force  transducer  (i.e.,  a strain 
ring) . The  subject  was  instructed  to  pull  on  the  rope 
using  only  his  neck  muscles,  to  build  rapidly  but  8SKX>thly 
to  a Euucimum  level,  and  to  hold  that  level  for  a count  of 
4 seconds.  The  subject's  feet  were  placed  flat  on  the 
floor  ana  the  subject  was  net  allowed  to  rise  up  from  the 
seat  or  use  his  torso  except  to  taaintain  his  position. 
Three  maximum  exertions  %#ere  made  in  each  of  the  four 
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directions  with  2 minutes  cf  rest  between  trials,  and  the 
average  fo*:ce  of  each  set  computed.  Figure  2.7  shows  typi 
cal  force  curves  and  the  EMG  signal  resulting  from  these 
tests. 


Figure  2.7  EMG  and  Force  Signals  Resulting  from 
Isometric  Strength  Tests 


2.  Results . T2ible  2,7  gives  the  strength  and  reflex 
time  statistical  results  for  the  18  Navy  sixbjects.  The 
abbreviated  variable  names  may  be  cross  referenced  in 
/‘.ppendix  A for  a more  complete  title. 


VArlABLE 


TABLE  2.7 

REFLEX  TIME  AND  STRENGTH  STATISTICS 
LE  H .1E1N  £U  DEV  aiMMUr.  ftAiitlUrt 


STF  E 
STF.  FLX.H 
S7FSAGAV 
iAGLATAV 


It  will  be  noted  that  the  average  reflex  times  for  the 
three  pull  directions  are  neaurly  the  seune,  being  between 
51  and  56  msec,  while  the  range  in  reflex  times  is  from  37 
to  70  msec. This  is  in  contrast  to  the  results  of  two  pre- 
vious studies  ( 22  and  23)  where  the  data  suggested  that 
the  reflex  tines  in  lateral  bending  were  significantly 
less  than  in  sagittal  bending.  The  strength  results  indi- 
cate nearly  identical  values  for  right  and  left  lateral 
pulls  and  these  are  nearly  the  seune  as  the  results  for 
flexion,  although  the  maximum  forces  achieved  in  lateral 
bending  are  significantly  higher  than  achieved  in  flexion 
(56.0  Ibf.  to  42.0  Ibf.).  The  greatest  strengths  were  in 
ejctension  where  the  average  of  46.0  Ibf.  is  about  33% 
greater  than  the  average  in  lateral  bend  or  flexion.  The 
maximum  strength  in  extension  was  similar  to  that  for 
lateral  bend,  however.  Table  B.4  in  Appendix  B shows  the 
individual  strength  and  reflex  values  from  which  these 
statistics  were  derived  where  each  strength  value  is  the 
average  of  three  trials  and  each  reflex  time  is  the  average 
of  at  least  6 tests. 

E.  Comparison  of  Measurement  Results  with  IIHS  Study 
Results' 

1.  Anthropometry . As  indicated  by  the  third  letter 
in  the  prefix  code  of  the  HSRI  subject  numbers  (S,M,  or  T) , 
eleven  of  the  eighteen  subjects  tested  are  of  medium 
stature,  six  are  tall,  and  only  one  is  short  according  to 
U.S.  population  data  on  stature  for  18-24  year  males. 

Thus,  the  group  of  subjects  used  in  this  study  is  biased 
toward  tall  individuals.  The  mean  stature  for  the  group 

is  177.0  cm  as  compar«d  to  174.86  cm  and  174.95  cm  for 
young  male  medium  stature  subject  groups  in  tha  sagittal 
and  lateral  studies  respectively. 

2.  Range  of  Motion.  Table  2.8  summarizes  the  average 
range  of  motion  results  for  the  primary  Euler  angles  in  the 
planar  head  movements  and  compares  these  results  with 


NIFARISON  OF  NAMRL  AND  TIHS  SIUDY  RANGE-0?’-?-iOTION  RESULTS 


average  results  from  the  different  subject  groups  in  the 
IIHS  lateral  plane  study.  It  can  be  seen  that  for  every 
position,  the  average  range  of  motion  for  the  NAMRL  sub- 
jects is  greater  than  for  all  other  groups  and  this 
difference  is  particularly  significant  in  extension.  As 
expected,  the  results  are  most  similar  to  the  18-24  year 
male  and  female  groups.  A t-te:<t  for  comparison  of  popu- 
lation means  shows  the  NAMRL  results  to  be  not  significant- 
ly different  from  the  18-24  year  male  results  at  the  .10 
level  of  significance.  It  is  intarerting,  however,  that 
the  standard  deviation  for  the  Euler  angles  is  consis- 
tently and  considcrcibly  smaller  for  the  NAMRL  population 
than  for  the  IIHS  study  population.  This  is  illustrated 
in  Table  2.9  which  compares  the  NAMRL  results  with  the 
results  from  18-24  year  miles  from  the  IIHS  lateral  study. 
Perhaps  this  is  due  to  the  difference  in  subject  motiva- 
tion, level  of  training,  and  subject  conditioning  which 
results  from  experience  as  subjects  for  huxoan  experiments. 

Table  2.9 

Comparison  of  Means  and  Standard  Deviations 
of  Range- of-Motion  Results  for  NAMRL  Subjects 
and  18-24  Year  Males  From  IIHS  Study. 

Primary  Euler  Angle  (degrees) 
IIHS  Lateral  Study  NAMRL 


Position 

Mean 

S.  D. 

Mean 

S.D. 

Extension 

72.8 

18.2 

79.0 

6.6 

Flexion 

56.2 

11.6 

6.  . 5 

5.4 

R.  Rotation 

73.2 

9.4 

76.1 

5.8 

L.  Rotation 

76.2 

7.2 

77.5 

5.9 

R.  Lateral  Bend 

41.7 

12.6 

47.1 

6.4 

L.  Lateral  Bend 

44.6 

11.3 

45.0 

7.5 
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3.  Reflex  Times  and  Muscle  Strenath.  Table  2.10 


compares  the  average  NAMRL  results  for  measured  muscle 
reflex  times  and  isometric  pull  force  with  .veraged 
results  from  the  IIHS  lateral  and  sagittal  studies.  With 
regard  to  reflex  times,  there  is  no  particular  group  that 
the  NAHRL  subjects  match  particularly  well  or  consistently 
for  all  planes.  As  mentioned  previously  the  NAMRL  data 
show  similar  reflex  times  for  all  directions  while  IIHS 
study  results  show  longer  times  for  sagittal  movements 
than  for  la*-eral  movements.  This  discrepancy  is  unex- 
plained at  this  time.  In  lateral  bending,  the  NAMRL 
results  match  best  with  the  62-74  year  females , the  18-24 
year  males,  and  the  35-44  year  males.  In  sagittal  bending, 
NAMRL  average  reflex  times  are  closest  to  the  young  and 
middle  aged  females. 

Concerning  strength  results,  it  is  seen  that  the 
NAMRL  subjects  match  extremely  well  with  the  35-44  year 
males.  A t-test  for  comparison  of  population  means  for 
these  two  groups  shows  no  significant  difference  at  a 
significance  level  of  .05. 
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Chapter  3 

COMPUTER  SIMULATIONS  OF  NAMRL  SLED  TESTS 

This  chapter  ia  divided  into  three  major  sections. 

In  Section  A,  a brief  description  of  the  characteristics 
and  capabilities  of  the  MVMA-2D  model  is  given.  In 
section  B,  these  characteristics  as  they  apply  to  this 
study  are  developed  fxirther  in  a description  of  how  the 
measurement  data  described  in  Chapter  2 have  been  used  to 
determine  model  par2Lineter  values  and  how  other  parameters 
for  which  there  are  no  data  available  vwre  determined. 
Section  C presents  the  simulation  results  for  6 and  15  G 
sled  tests  and  is  divided  into  two  parts.  The  first  part 
describes  simulation  results  obtained  by  using  experimental 
Tj^  acceleration  signals  as  direct  input  to  T^.  In  the 
second  part,  results  obtained  by  using  sled  acceleration 
data  and  including  restraint  system  and  torso  characteris- 
tics are  presented. 

A.  The  MVMA-2D  Crash  Victim  Simulator 

All  simulation  work  in  this  study  has  utilized  the 
MVMA-2D  Crash  Victim  Simulator,  Version  III.  This  model, 
in  its  current  form,  is  a result  of  H5RI  extensions  and 
improvements  made  upon  the  original  CAL  2-D  model  (1966) , 
later  modified  to  the  ROS  (Revised  Occupant  Simulation)  in 
1971  and  MODROS  (Modified  Revised  Occupant  Simulation)  in 
1972.  Pertinent  to  the  modeling  in  this  study  is  the  use 
of  a two- joint  extensible  neck.  The  human  occupant  is 
constructed  of  9 body  segments  (head,  neck,  upper  torso, 
middle  torso,  pelvis  or  lower  torso,  upper  log,  lower  leg, 
upper  arm,  and  lower  arm)  divid*Ki  by  8 pivot  joints  as 
shown  in  Figure  3.1.  For  each  segment  except  the  neck, 
the  mass,  location  of  the  center  of  mass,  and  mass  moment 
of  inertia  are  specified.  For  the  neck,  the  mass  is  dis- 
tributed as  desired  at  the  upper  and  lower  neck  pivot 
points . 
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The  equation  of  motion  describing  the  dynamic  be- 
havior of  the  articulated  occupant  were  derived 
using  Lagrangian  formulations.  Energy  dissipation  at  the 
joints  may  be  through  the  mechanisms  of  friction  and/or 
viscosity  while  the  moments  acting  at  each  of  the  simulated 
joints  may  be  derived  from  up  to  five  sources  including 
biodynamic  muscle  tension,  elasticity,  viscous  damping, 
coulozob  friction,  and  non-linear  energy-dissipating 
motion-restrictive  stops.  Interaction  of  the  occupeuit 
with  vehicle  structures  and  restraint  systems  may  be  es- 
taiblished  by  specifying  contact  ellipses  with  desired 
material  properties  (including  force  deformation  and  energy 
^d}sorbing  characteristics)  and  contacting  surfaces  and 
belts  with  specified  properties.  The  restraint  system  can 
utilize  up  to  four  belts  — two  attaching  to  the  hip 
segment,  one  to  the  upper  torso,  and  one  attaching  arbi- 
trarily to  any  torso  segment.  In  addition,  options  for 
free  slip  between  belt  pairs  (force  equalization)  or  for 
percentage  force  limits  of  one  belt  relative  to  another 
for  simulating  friction  are  also  available.  A more  de- 
tailed description  of  these  features  may  be  found  in 
references  2 and  20. 

B.  Determination  of  Model  Parameters 

A valuable  feature  of  the  KVMA-2D  model  as  it  relates 
to  this  study  is  that  it  has  been  developed  based  upon 
attempts  to  consider  and  simulate  the  individual  physical 
factors  which  have  an  influence  on  the  dynamic  response 
of  the  occupemt.  For  exeimple,  instead  of  lumping  all 
neclc  properties  into  a general  visco-elastic  element  with 
two  parameter  constants,  the  model  is  general  enough  to 
allow  separation  of  some  of  the  factors,  such  as  muscle 
versus  passive  tissue,  which  physically  act  to  affect 
the  response.  Thus,  it  is  ix>ssible  to  use  experimental 
data  where  available  to  simulate  these  factors  and  to 
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gain  a better  understanding  through  the  model  of  the  re- 
lative importance  of  each. 

At  the  beginning  of  the  modeling  phase  of  the  project, 
it  was  decided  to  match  simulation  results  with  the  sled 
test  results  from  a closely  matched  subgroup  of  the  18  sub- 
jects rather  than  to  match  with  all  18  subjects  or  with  an 
individual  subject.  An  examination  of  the  emthropometri c 
and  strength  measurements  for  those  properties  judged  to 
have  the  most  significant  effects  on  head  and  neck  dynamic 
response  in  sagittal  flexion  was  made  and  resulted  in  a 
aelectio.i  of  5 subjects  (NAMOl,  NAM04 , NAH06,  NAM08,  and 
NATlS)  whose  measurement  data  wouid  be  used  for  establish- 
ing model  parauneter  values  and  whose  sled  test  results 
would  be  used  for  comparison.  Tables  3.1  and  3.2  give  the 
individual  data  and  group  statistics  for  selected  measure- 
ments used  in  establishing  the  data  set  of  the  NAMRL  simu- 
lations. While  these  results  were  uoed  to  establish  model 
parameter  values,  tb“  manner  in  which  these  data  are  used 
depends  upon  certain  modeling  assumptions  and  other  non- 
obvious  procedures.  In  addition,  for  other  parameters 
used  in  the  model  there  are  no  experimental  data  from 
which  to  derive  reliable  parameter  constants.  For  these 
reasons,  the  following  subsections  are  included  to  docu- 
ment the  procedures  used  to  date  in  this  study. 

1 . Segment  specifications  — lengths,  masses, 
centers  of  masses,  and  moments  of  inertia. 

a.  Torso  ana  Extremities.  Tox'so  length  was 
computed  as  the  distance  from  trocanterion  height  to  cer- 
vicale  height  as  measured  in  the  seated  position.  As 
shown  in  Figure  3.1,  the  torso  is  divided  into  three 
segments.  Initially,  the  lengths  of  these  individual 
segments  were  de+'ermined  by  proportioning  them  the  sesme 
relative  to  total  torso  length  as  is  used  in  the  MVMA-2D 
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TABLE  3,1  SELECTED  ICEASUREMEKTS  FOR  FIVE  KAMRL  SUBJECTS 


TABLE  3.2 

S£L£CTIt  KUSOBEISFJII  oTA^ISTlCS  FOi-  5 « AVT  SiJxiJECTS 


¥ AEIABIF 

K 

91AK 

SIC  CBV 

KIBIAU9 

BAXiaua 

BT  (KG) 

5 

73. C 

4.  1 

65.7 

75,5 

STAT  (CM) 

5 

176.9 

2.9 

174.2 

181.  7 

E9SITBT 

5 

92. S 

1.  3 

91 . 6 

94.  9 

tZATCIF 

5 

57.7 

2.7 

54.7 

61.  6 

FACiH-: 

5 

13.3 

C.  a 

11.9 

14.  1 

AEHK3B 

5 

11.0 

0.5 

10.4 

11.  7 

SUFMKCIE 

3 

37.2 

1.2 

35.4 

38.7 

lUPliKClF 

5 

36.0 

1.  1 

36.5 

39.  3 

ACEEADLG 

5 

32.1 

1.3 

3(i . 5 

34.  1 

SADSTTIG 

5 

25.  a 

0.7 

24.9 

26.  7 

HANDLG 

5 

18.5 

C . 6 

18.  J 

19.  9 

TP.Cf  2H1G 

5 

40.7 

1.4 

39.3 

4 2.  8 

f IBOLAIG 

5 

40.7 

1.5 

39.5 

43.  2 

liPESITHT 

5 

89.5 

1.5 

87.7 

91.  4 

tbagb:  5 

c; 

75.6 

1.4 

73.9 

77.  3 

GLABLBIS 

5 

73.7 

1.6 

78.3 

82.  4 

C7HT  S 

5 

64.6 

1.7 

63.1 

67.  3 

Tsocnais 

5 

9.6 

0.3 

6.9 

10.  9 

C7  SX 

5 

-9.2 

2.6 

-12.4 

-7.  1 

C7  SZ 

5 

€4.6 

1.7 

€3.  1 

67.  3 

T9AG  SX 

5 

-1.2 

3.4 

-4.V 

2.  5 

TFAG  SZ 

c 

75.6 

1.4 

73.^ 

77.  3 

BFL  A¥G 

5 

56.5 

6. 3 

50.  3 

63.  3 

STB  2XTF 

5 

46.3 

7.6 

37.  J 

55.  0 

STB  FLXB 

5 

36. C 

2.3 

33.  3 

39.  3 

F2KE0IE 

t 

0.6 

2.5 

-3.1 

3.7 

F3EXIE 

5 

83.0 

4.2 

76. 3 

87.  5 

BUFLFXE 

3 

-61.5 

7.3 

-68.2 

-50.  6 

FSAGSCK 

5 

144.5 

9.0 

134. 3 

155.  7 
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baseline  data  set^  for  5Cth  percentile  males.  In  this  data 
set  the  lower-most  segment  corresponds  to  the  pelvic  mass 
and  is  20%  of  the  torso  length.  The  middle  and  upper  torso 
segments  are  about  25  and  55  percent  of  the  torso  length 
respectively.  These  proportions  are  somewhat  arbitrary,  how- 
ever, and  were  later  changed  in  this  study  to  44  and  36  per- 
cent respectively  in  order  to  provide  for  torso  bending 
above  the  chest  restraint  belt.  In  most  of  the  results 
presented  later  in  this  chapter,  however,  the  upper  and 
middle  torso  joints  (joints  3 and  4)  have  been  made  essen- 
tially rigid  so  that  the  torse  bends  only  at  the  hip  joint 
(see  section  7) . 

Extreinity  segment  lengths  were  determined  directly  from 
the  averaged  anthropometric  measures  taken  in  this  study. 
Initially,  trocanter-f emoral  length  was  used  for  upper  leg 
length,  fibula  length  for  lower  leg,  radiale-stylion  plus 
hand  length  for  lower  arm,  and  acromion-radiale  for  upper 
arm  length.  At  a later  point  in  the  study,  the  arm  segments 
were  removed  and  one  half  of  the  upper  arm  :r>ass  was  added  to 
the  upper  torso  mass.  This  was  done  when  it  was  realized  from 
the  high  speed  films  that  the  arms  were  restrained  by  straps. 

Mass  and  moment  of  inertia  values  for  the  torso  were 
also  scaled  in  proportion  to  segment  values  in  the  baseline 
data  relative  to  total  body  mass  minus  the  mass  of  the  head 
and  neck.  Segment  masses  were  scaled  in  direct  proportion 
to  segment  length  while  moments  of  inertia  were  scaled  to 
the  oaseline  data,  by  proportions  of  mass  times  segment 
length  squared. 

Distances  fro.m  link  ends  (i.e.,  joints)  to  segment 
centers  of  mass  were  scaled  to  the  baseline  data  propor- 
tions. These  values  were  all  adjusted  appropriately  when 


This  baseline  date  set  for  the  MVMA-2D  model  has  been  es- 
tablished and  i.iodified  over  the  years  from  existing  and  newly 
acquired  data  and  continues  to  be  updated  and  improved  as  the 
mexial  is  used. 
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the  upper  torso  joint  (joint  3)  was  moved  up  . The 
values  shown  in  Table  3.3  are  those  used  in  the  NAMRL 
data  set  for  the  results  of  this  report. 

Table  3.3 


Torso  and  Extremity  Segment  Specifications 
For  NAMRL  Data  Set 

End  of  Link  to 

Segment  Length (cm)  Canter  of  Mass (cm)  Mass (kg) 

iyy  2 

(Kg-m  ) 

Upper  torso 
(includes 
1/2  upper 
arm  mass) 

20. 3 

10.15 

15.3 

.1433 

Middle  torso 

23.9 

8.8 

9.5 

.1343 

Hip 

10.8 

4.2 

8.4 

.1995 

Upper  leg 

40.7 

18.5 

17.62 

.2705 

Lower  leg 

40.7 

28.8 

9.5 

. 3412 

b.  Head  and  Neck  Mass  and  Moment  of  Inertia.  h 

correlation  of  anthropometric  dimensions  with  head  mass 

aud  moment  of  inertia  measurements  on  five  m.ale  cadavers 

from  a study  by  Chandler,  et  al.,  1975  (3)  showed  that 

head  loass  is  highly  correlated  with  head  circumference  and 

that  moment  of  inertia  is  highly  correlated  with  the 

2 2 

quantity  [ (menton  to  vertex)  + (head  length)  ] x [head 
circumference].  Accordingly,  these  anthropometric 
measures  from  the  five  NAMRL  subjects  were  used  with  these 
measures  from  cadavers  to  obtain  estimates  for  head  mass 
and  head  moment  of  inertia.  Neck  mass  was  obtained  by 
scaling  to  the  head  mass  in  proportion  to  tae  MVMA  base- 
line daca,  and  was  distributed  with  33  percent  at  the 
condyles  and  67  percent  at  the  lower  neck  joint.  To  the 
head  mass  and  moment  of  inertia  were  added  the  mass  and 

'omcnts  of  inertia  due  to  the  instrument  packages  riven 

2 

as  .522  kg  and  .0075  kg-m  respectively  by  Ewing  and 
Thomas  (9  ).  These  calculations  resulted  in  the  following 


42 


values  for  the  NAMRL  data  set. 


Instrumented  Head  Mass  = 4.615  Kg 

2 

Instrumented  Head  lyy  = .024  Kg-m 

Neck  Mass  = 1.194  Kg 

c.  Neck  length  and  location  of  head  c.g.  While 
x-rays  of  the  NAMRL  subjects  were  not  available,  an 
estimate  of  neck  length  and  location  of  the  head  center  of 
gravity  relative  to  the  occipital  condyles  was  obtained 
by  using  x-ray  films  from  young  male  subjects  of  the  IIHS 
studies.  The  X and  Z distances  from  tragion  to  the  con- 
dyles were  measured  and  scaled  and  added  to  the  average 
distances  of  tragion  to  head  center  of  gravity  determined 
by  Ewing,  et  al.  ( 8)  (i.e.,  the  head  c.g.  lies  2.1  cm 

above  and  1.3  cm  foirward  of  tragion).  To  this  was  added 
the  distance  the  head  c.g.  is  shifted  by  the  instnomen- 
tation  package  which  is  given  by  Ewing  and  Thomas  (9 ) as 
.35  err.  forvard  and  ,2  cm  down.  Neck  length  was  estimated 
by  using  anthropometry  results  for  the  five  NAMRL  subjects 
in  order  to  locate  tragion  and  cervicale  in  two  dimensions. 
X-ray  measurements  from  the  IIHS  young  males  were  utilized 
in  order  to  locate  tlie  condyles  and  C^-T^^  relative  to 
these  external  anatomical  points.  The  neck  length  was  then 
computed  as  the  straight  line  distance  from  C^-T^  to  tlie 
occipital  condyles.  The  results  of  these  measurements  and 
calculations  gave  the  following  values: 

X distance  condyles  to  head  c.g.  = 2.47  cm 
z distance  condyles  to  head  c.g.  =-4.16  cm 
neck  length  =11.2  cm 

where  the  positive  x axis  is  forward  in  the  Frankfort 
plane  and  positive  z is  down. 


2.  Head  and  Neck  Range  of  Motion.  Since  the  physical 
situation  being  simulated  involved  priri\arily  sagittal  plame 
movements^  only  the  sagittal  plane  range  of  motion  results 
need  to  be  considered  here.  These  results  give  a measure 
for  the  maximum  head  aiigle  forward  (flexion)  and  maximum 
head  angle  rearward  (extension)  from  the  Frankfort  plane 
position  achieved  by  voluntary  effort.  While  these  angles 
are  the  cumulative  result  of  bending  at  several  articula- 
tions along  the  length  of  the  neck  (i.e.,  at  each  cervical 
disk)  and  at  the  condyles,  the  MVMA-2D  model  considers  only 
two  neck  joints  connected  by  a straight-line  segment  neck. 
While  these  joints  may  be  positioned  as  desired,  it  was 
considered  most  reasonable  to  initially  consider  the  upper 
neck  joint  to  be  at  the  occipital  condyles  and  the  lower 
neck  joint  to  be  at  the  C^-T«  disk.  For  range-of-motion 
input  specif icationc , then,  the  model  requires  that  joint 
stop  angles  for  movement  of  the  neck  relative  to  the  torso 
and  movement  of  the  head  relative  to  the  neck  be  specified. 
Figure  3.2  illustrates  these  required  stop  angles  where  the 


Figure  3.2  Ranne-of-Motion 
ZD  Crash  Victim  Simulatcr. 
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maximum  flexion  and  extension  of  the  neck  relative  to  the 
torso  are  le±>eled  y and  6 respectively,  and  the  maximum 
angle  of  the  head  beyond  the  neck  angle  in  full  flexion 
and  extension  are  a and  6 respectively.  The  problem,  then, 
is  to  determine  these  angles  from  the  range-of-nvotion 
results.  Two  possible  solutions  were  considered. 

The  first  approach  considered  was  to  measure  a,  g,  y, 
and  6 directly  from  y-camera  photographs.  The  difficulty, 
of  course,  is  to  locate  the  condyles  and  the  disk. 

These  points  were  estimated,  however,  by  using  ^iverage 
measures  of  the  x and  2 distances  from  C^-T^  to  cervicale 
emd  tragion  to  the  condyles  obtained  from  x-rays  of  young 
males  in  the  IIHS  study.  These  distances  were  then  scaled 
appropriately  and  marked  off  from  cervicale  and  tragion  on 
the  projected  NAMRL  range-of -motion  photographs.  In  this 
way,  the  angle  of  the  neck  relative  to  the  torso  lire 
(assumed  to  be  vertical)  was  measured  directly  in  flexion, 
extension,  and  Frankfort  positions.  The  angle  of  the  head 
relative  to  the  neck  was  measured  as  the  angle  between 
the  perpendicular  to  the  Frankfort  plane  and  this  neck 
angle  line.  While  there  is  some  inaccuracy  in  the 
measurements  due  to  the  difficulty  in  locating  cervicale 
in  extension  and  the  uncertainty  of  t}ie  change  in  orien- 
tation of  Cy  during  flexion  and  extension,  some  consis- 
tency in  the  method  was  found  and  average  values  for  a, 
g,  Y,  and  i were  determined  to  be: 

u = -3.5“ 

g - 64.9° 

Y = 64.4° 

6 = 19.0° 

The  sum  of  these  angles  should  equal  the  total  sagittal 
range  of  notion  and  is  144.8  degrees  compared  to  144.5 
degrees  determined  by  orthogonal  photogrammetry  techniques 
for  these  five  subjects. 


When  the  joint  stop  angles  are  determined  in  this 
way  (where  the  orientation  of  the  neck  is  considered  to 
be  that  of  the  straight  line  segment  connecting  to 

the  condyles)  some  of  the  ramge  of  motion  due  to  articu- 
lations at  neck  joints  is  attributed  to  the  upper  neck 
joint  in  the  model.  While  this  may  be  a reasonable  way 
of  dividing  the  range  of  motion  between  the  two  joints, 
it  does  have  some  drawbacks.  If  the  upper  neck  joint  is 
to  represent  the  occipital  condyles,  which  is  an  important 
articulation  between  the  large  mass  of  the  head  and  the 
relatively  small  extensible  neck,  then  it  is  perhaps  more 
important  that  the  joint  stop  angle  specification  for  this 
joint  be  correct  even  though  this  means  that  the  lower 
joint  assumes  tha  cumulative  range  of  motion  for  all  other 
neck  joints. 

As  a result  of  these  considerations,  a second 
approach  was  used  to  determine  the  joint  stop  angles  needed 
in  the  model.  In  the  IIHS  sagittal  study,  x-rays  were 
taken  while  subjects  performed  full  flexion  and  full  ex- 
tension movements.  From  these  x-rays,  the  change  in 
angle  of  C2  (this  is  nearly  the  same  as  the  angle  of 
and  easier  to  measure)  relative  to  the  vertical  was  used 
for  the  lower  neck  joint  range  of  motion,  while  the  change  in 
head  angle  (as  determined  by  the  Frankfort  plane)  relative 
to  was  used  for  the  condyle  range  of  motion.  The 
results  of  these  measures  gave: 

a = -5.0° 

a - 25.5° 

Y = 71. C° 

6 = 53.3° 

While  these  values  were  not  determined  from  the  NAMRL 
subject  data,  the  total  range  of  motion  which  is  the  sum 
of  these  angles  is  144.3  degrees  compared  to  144.5  for  the 
fi’  a NAMRL  subjects.  Also,  a t-test  between  the  IIHS 
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sagittal  study  range-of-motion  results  and  the  NAMRL  range- 
of-motion  results  showed  no  significant  difference  at  the 
.10  level  of  significemce. 

The  primary  difference  in  these  angles  from  those 
derived  by  the  first  method  is  in  e^ctension  where  the 
head-neck  angle  is  25.5**  compared  to  64.9°  and  the  neck- 
torso  angle  is  53.3“  compared  to  19“ . For  the  first 
methody  the  total  head  flexion  from  vertical  is  a -f  r ~ 60.9“ 
while  for  the  second  it  is  66“ . In  extension^ the  first 
method  gives  e + 6 » 83.9“  as  ccmpared  to  76.3“  for  the 
second  method. 

While  NAMRL  simulations  using  stop  angles  determined 
by  both  of  these  techniques  have  been  used,  simulations 
presented  in  this  report  use  the  results  of  the  second 
approach.  Differences  in  occupant  response  due  to  the  two 
sets  of  data  depend  upon  stiffness  values  of  the  joint 
stops,  however.  As  will  be  discussed  in  the  next  section, 
the  values  for  joint-stop  stiffness  being  used  at  this 
time  result  in  only  small  differences  in  the  model  output 
from  the  two  sets  of  data. 

• Passive  Joint  Torques  and  Joint  Stops.  There 

are  virtually  no  data  available  from  which  to  determine 

reasonable  estimates  for  passive  joint  resistance  and 

joint  stop  characteristics  of  the  head  and  neck.  While 

the  MVMA-2D  model  has  the  eapahility  to  simulate  these 

torques  by  several  sources,  initially  all  jjassive 

resistance  coefficients  within  the  range-of-motion  angles 

were  set  to  zero  and  torques  due  to  movement  beyond  the 

joint  stops  were  represented  by  a constant  times  the 

square  of  the  angular  deformation  of  the  stop.  Values 

which  have  been  used  previously  in  the  MVMA-2D  model  for 

this  quadratic  angular  deflection  coefficient  are  2.82 

2 

and  5.5  N-iVdeg  for  both  neck  joints  at  flexion  and 
extension  stops  respectively.  When  the  model  was  run 
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under  these  conditions  (with  muscle  effects  included) , 
however , the  results  showed  sharp  peaks  of  acceleration 
producetd  by  tlie  joint  stops  which  were  clearly  not  present 
in  the  experimental  sled  test  results  (see  Figure  3.23). 

It  was  clear  that  the  joint  stops  were  too  stiff  but 
intuitively  the  values  used  did  not  seen  excessive  for  the 
range-of -motion  angles  measured  on  the  subjects. 

An  attempt  to  remedy  this  problem  was  made  by  reducing 
the  joint  stop  angular  deflection  coefficients  by  a factor 
of  100  and  also  adding  some  apprecizd:)le  viscous  damping  to 
the  joints  throughout  the  range  of  motion.  The  results 
were  much  more  reasonable.  While  little  justification  can 
be  offered  at  this  time  for  these  values  other  than  the 
fact  that  good  results  are  obtained,  the  expleination  per- 
haps lies  in  the  inappropriatenesF  of  the  term  "joint 
stops,"  For  the  purposes  of  this  study,  the  joint  stop 
angles  were  defined  to  be  those  angles  at  which  a human 
subject  is  unable  or  unwilling  to  move  his  head  any 
further  in  a given  direction.  No  measurements  have  been 
made,  however,  of  the  forces  it  took  to  achieve  this  final 
position.  There  is  little  doubt  that  for  motivated  sub- 
jects the  torques  required  to  go  beyond  their  voluntary 
effort  are  compatible  with  tlie  magnitude  of  the  stiffness 
coefficients  used  in  the  MVMA-2D  model,  and  that  the 
assumption  that  the  torque  is  proportional  to  the  square  of 
the  angle  beyond  the  stop  is  reasonable.  But  it  is  also 
reasonable  to  expect  that  the  joint  torques  have  been 
increasing  in  some  manner  up  to  this  final  position.  By 
reducing  the  quadratic  coefficients  and  adding  the  viscous 
damping  constant,  it  was  possible  in  some  suitable  but 
still  inexact  way,  to  model  the  complex  characteristics  of 
the  physiological  joint.  In  addition,  by  adding  the  vis- 
cous damping  coefficient,  the  resistive  torque  is  made 
sensitive  to  angular  velocity  which  also  seems  to  be  a 
reasonable  attempt  to  model  effects  of  passive  tissue  such 
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as  ligaments.  There  is  no  doubt,  however,  that  further 
research  and  study  is  needed  in  this  area  before  confidence 
in  this  aspect  of  the  model  will  be  achieved.  In  any  case, 
for  the  results  presented,  neck  viscous  damping  coeffi- 
cients have  been  set  to  .01  and  .03  N-m-sec/deg.  for  the 
upper  and  lower  neck  joints  respectively.  For  the  joint 
stop  quadratic  deflection  coefficients  which  become 
effective  when  the  joint  stop  angles  are  exceeded,  the 
following  values  have  been  used: 

K (flexion  and  extension)  = .0087  N-m/deg^ 

^condyles  . .0261  N-m/deg^ 

^condyles  (extension)  = 1.0  N-m/dog^ 

In  choosing  the  value  for  the  condyles  in  extension,  con- 
sideration was  given  to  the  fact  that  the  subject  is 
initially  positioned  very  close  to  his  measured  head-neck 
joint  stop  angle  in  extension  isee  sections  2 and  6)  where 
the  resistive  torques  are  significant.  Consequently,  this 
constant  was  maintained  reasonably  close  to  the  MVMA 
baseline  value. 

4.  Neck  Muscle  Model.  There  arc  several  models  in 
the  literature  which  attempt  to  simulate  muscle  based  upon 
experimental  observations,  all  of  which  have  limitations 
and  deficiencies  and  are  simplified  approximations  of  a 
complex  mechanism.  The  MVMA-2D  simulator  models  the  active 
element  by  a spring  and  dashpot  in  series  as  shown  in 
Figure  3.3. 

c(M) 

k(M)  n 

• AA/V 1 


Figure  3.3  Muscle  Element 


Figure  3,4  Simplified  Free-Body  Diagram  of  Head  and 
Neck  Showing  Major  Forces  Involved  During  Isometric 
Strength  Testing . 
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Here  the  coefficients  K(M)  and  C{M)  are  considered  simple 
functions  of  the  voluntary  static  moment,  M (i.e.,  of  the 
"tightness"  of  the  muscles),  as  given  by  the  following 
equations : 

K ■ ai  + ao  |Mj 

C = aajHi 

From  these  equations  it  is  seen  that  when  the  muscles  are 
relaxed  (i.e.,  M » 0) , the  murcle  has  no  effect  on  joint 
torque  since  C - 0.  The  muscle  tension  is  time-dependent 
and  is  input  to  the  model  in  tabular  form. 

The  values  of  aj,  a2,  and  33  in  the  above  equations 
are  joint  parameters  and  are  dependent  on  the  particular 
muscle  strengths  of  the  occupant  and  the  particular  joints 
involved.  Baseline  data  for  these  values  have  been  deter- 
mined by  Bowman  (1)  and  are  derived  from  experimental  data 
on  the  knee  joint  obtained  by  Moffatt,  Hassis,  and  Haslam  (17) . 
For  the  lower  neck  joxnt  and  |H{  max  * 27.73  N-m,  Bo%raaan  gives; 

aj  ■ .1668  N-m/deg 
- 1 

32  = .153  deg 

33  “ .0129  sec/deg 

For  the  NAMRL  population,  the  actual  forces  in  the  neck 
muscles  at  the  upper  and  lower  neck  joints  were  estimated  from 
the  measured  isometric  pull  forces  by  summing  moments  about 
the  condyles  and  Cy-T,  respectively  as  shown  in  Figure  3.4. 

The  distance  £2  estimated  from  measurements  taken  from 

tragion  to  the  head  band  during  strength  tests  and  by  x-ray 
measurements  from  the  condyles  to  tragion.  Average  values 
used  for  the  NAMRL  population  were  3.32  cm  and  2.1  cm  res- 
pectively giving  an  £2  distance  of  5.42  cm.  The  value  for 
was  the  neck  length  as  computed  in  Section  1.  For  the 
extensor  musclas,  £;  and  £3  were  estimated  at  2"  or  5.08  cm. 
from  x--ray  and  skull  measurements. 
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Using  these  values  and  the  average  muscle  pull  force 
for  the  five  subjects  in  extension  (48.28  Ibf.),  the  values 
of  T_  and  T'  were  calculated  as: 

» 51.5  Ibf.  » 229.1  N 

£ 

T'  = 157.96  Ibf.  * 702.6  N 


While  the  baseline  values  for  a^,  a2 , and  as  were 
computed  assuming  a maximum  joint  torque  of  27.73  N-m  due 
to  muscle,  these  values  were  used  for  the  lower  neck  joint 
muscle  model  where  the  maximum  estimated  torque  in  extension 
is  5.08  cm  times  702.6  N or  35.7  N-m  and  they  were  scaled  to 
the  condyle  joint  where  the  maximum  estimated  torque  in 
extension  is  11.6  N-m.  Bowman  (1)  has  shotim  that  a^,  a^ , and 
a3  may  be  reasonably  sealed  as  follows  from  joint  to  joint  in 
an  individual  or  from  individual  to  individual: 


where  th<’  subscrifJts  1 <in<J  II  refer  to  the  two  joints  involved. 
Since  the  distances  through  which  thr  extensors  act  are  assumed 
the  same  for  both  neck  joints,  the  only  constant  that  is 
altered  for  the  condyles  is  aj.  This  is  scaled  by  the  ratio 
of  the  maximum  muscle  tensions  which  is: 


T 


E 


2 3C.0  N 
702.6  N 


- 326 
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Table  3.4  gives  the  muscle  parameter  constants  used  in 
simulations  of  this  report  based  upon  these  assumptions  and 
calculations. 


Table  3.4 

Muscle  Pareuneter  Values  Used  in  NAMRL  Simulations 


Parameter  

10b%  Muscle 

Joint  a;(N-m/deg)  a;^.  (deg  ')  83  (sec/deg)  Torque(N-m) 

Occ.  Condyles  .053  .153  .0129  11.6 

C7-T1  ,1668  .153  .0129  35.7 

For  neck  stretch,  these  constants  were  converted  to 
lineal  coefficients  (as  opposed  to  angular)  and  scaled  appro- 
priately by  considering  the  total  muscle  tension  due  to  flex- 
ors and  extensors.  For  computing  flexor  muscle  tension  a 
similar  technique  of  taking  moments  about  Cy-T;  was  used  con- 
sidering the  primary  muscle  group  (the  sternomastoids)  to  act 
at  an  average  distance  of  about  1"  (i^)  relative  to  the  lower 
neck  joint.  For  the  average  flexion  pull  strength  for  the  5 
NAMPL  subjects  of  36  Ibf.  and  the  same  values  of  and  as 
for  extension  (see  Figure  3.4),  a maximum  flexor  tension  of 
235.5  Ibf,  or  1047.5  N is  calculated.  When  added  to  the  maxi- 
mum extensor  muscle  tension  of  702.6  N,  the  total  neck  muscle 
tension  is  estimated  to  be  1750.1  N or  393  Ibf.  Bowman  (1) 

has  shown  that  the  muscle  parameter  constants  in  neck  stretch, 
£ £ t 

ai,  ay*  and  a3,  are  related  approximately  to  the  angular 
coefficients  used  above  — a;,  ay,  a-.  — by: 


5 3 


lihiTTk 


(4) 


® 2 * 


2A2/1 


I 2 

ai  a 2ai/A  ; 

af  “ a3/£ 


where  I is  the  moment  arm  distance  for  the  angular  constants « 
taken  in  this  study  to  be  the  average  for  extension  and 
flexion  or  1.5".  Using  a],  a^#  and  33  from  the  baseline  data 
gives : 


£ 

* 131.7 

N/cm 

1 

A2 

= 2.3 

- 1 
cm 

1 

a3 

* .193 

sec/cm 

In  addition  to  these  parameter  values  which  the  model 
uses  to  compute  K(m)  and  C(m),  values  for  the  maximum  muscle 
tension  and  torque  are  required.  In  previous  studies  (22,23) 
it  has  been  assumed  that  for  the  situation  where  there  is 
prior  warning  of  impact,  an  individual  will  be  able  to  pre- 
tense his  neck  muscles  to  100%  or  more  of  the  tension 
measured  in  an  isometric  laboratory  exertion.  On  further 
examination  of  this  question  in  this  study,  it  appears  that 
this  may  be  an  inaccurate  assumption.  Based  u{>on  subiective 
feelings  <incl  a brief  si  miy  <>l  KMC  ::i<|ital.s,  it  apjM'ars  that  an 
individual  is  only  able  to  develop  about  1/4  to  1/2  of  his 
maximum  isometric  pull  strength  when  tensing  without  an  ex- 
ternal reacting  surface.  Thus,  even  for  the  NAMRL  subjects 
who  are  fully  prepared  fur  the  impact  and  are  aware  of  the 
exact  time  it  will  occur,  it  is  questionable  whether  they 
can  be  "fully"  twsed  at  time  t = 0.  As  a result  of  these 
considerations,  muscle  tensions  and  torques  were  set  at  33% 
of  maximum  and  maintained  constant  throughout  the  simulations. 
Whether  the  constant  tension  assumption  is  reasonable  is  un- 
known at  this  time.  It  may  be  that  the  muscles  build  quickly 
(in  50  msec  or  less)  to  their  maximum  tension  once  the  head 
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begins  to  angulate,  and  this  time  may  be  even  shorter  than 
the  reflex  times  measured  in  this  study  since  the  muscles  are 
already  in  an  active  state.  It  may  also  be,  however,  that 
the  muscle  is  somehow  "cut  out"  by  the  violent  stretching  pro- 
duced by  the  impact.  The  problem  requires  further  study  and 
experimental  testing. 

In  any  case,  rather  than  attempting  to  guess  some  time- 
dependent  muscle  input  function,  a constant  muscle  tension  or 
torque  of  33%  of  maximum  was  used.  This  resulted  in  the 
following  values  for  the  five  NAMRL  subjects: 

Maximum  muscle  torque  about  condyles  = 3.87  N-m 

Maximum  muscle  torque  about  C7-T1  = 11.9  N-m 

Maximum  muscle  force  in  neck  stretch  = 583.4  N 

5 ♦ Neck  Stretch  and  Compression  imeters.  The  MVMA- 
20  model  simulates  the  stretch  and  compression  characteristics 
of  the  neck  due  to  passive  tissue  by  a spring-damper  combination, 
von  Gierke  (24)  has  reported  that  the  undamped  natural  frequency 
of  the  head  caused  by  z-excitation  of  the  upper  torso  is  about 
30  Hertz.  A first-order  spring  rate  for  the  neck  can  be 
approximated  by: 

= 4ir'f  ^ (M.  + 1/3  M)  ; 

s u n n 

where  = head  mass  and  M^  = neck  mass. 

For  NAMRL  subjects  Mj^  =*  4.093  Kg  and  = 1.533  Kg.  Using 
f^  = 30  Hz  gives  * 1636  N/cm. 

For  a mass-spring-damper  model  of  the  human  body  with 
spinal  column,  von  Gierke  (24)  gives  a range  of  .221  to  .266 
for  the  damping  ratio  for  the  composite  spinal  column.  For 
lack  of  better  data,  an  estimate  of  one  fourth  of  this  has  been 
assumed  reasonable  for  the  cervical  spine  alone. 
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Using  a value  of  .243  from  von  Gierke’s  data,  and  with  the 
critical  damping  value  given  by  2/  + 1/3  M^)  the 

damping  coefficient  is  given  by: 

C 

l/4  (.243)  = — - ^ — 

2/  K(M.  + 1/3  M„) 
n n 

C = .0608  X 2/1636(4.093  + .511) 

s 

C =10.55  N-sec/cm 

s 

6.  Neck  and  Head  Initial  Angles.  Initial  neck  angle 
was  based  upon  the  results  of  Ewing,  et  al.,  1975  (10)  in  which 
the  neck  angle  was  calculated  from  the  coordinates  of  the  head 
anatomical  and  ■'.’’i  anatomical  origin  locations  at  first  sled 
motion.  By  this  procedure,  the  neck  angle  is  estimated  by  a 
line  drawn  from  the  anterior- superic  ' corner  of  Tj  through 
tragion.  For  the  five  NAMRL  su'/jects  the  average  neck  angle 

is  approximately  20°  to  the  vertical. 

Head  angle  is  determined  by  the  pitch  orientation  of  the 
Frankfort  plane  relative  to  the  vertical.  Initial  head  angle 
was  determined  from  the  sled  test  experimental  curves  for  head 
pitch  angle  at  time  t = 0.  For  the  6 G runs  the  average  head 
angle  for  the  five  subjects  was  calculated  to  be  95.54  degrees 
(head  back)  while  for  the  15  G runs  the  average  head  angle  was 
calculated  as  93.64  degrees  (head  back). 

7.  Restraint  System.  One  of  the  most  difficult 
aspects  of  the  sled  test  simulation  was  Ln  modeling  the 
restraint  system.  While  complete  satisfaction  in  this 
area  has  still  not  been  achieved,  some  reasonable  progress 
has  been  made  as  reported  in  the  results  in  Section  C of 
this  chapter.  As  an  initial  step,  two  pieces  of  the 
webbing  niaterial  used  in  the  restraint  system  were  tested 
for  stress-strain  characteristics  under  static  loading 


conditions.  From  these  tests  a value  of  K « 112,000  N 
per  unit  strain  was  determined  and  used  in  the  model. 

Because  of  the  manner  in  which  the  belts  are  placed 
during  the  sled  tests  and  because  actual  belt  loading 
curves  are  not  measured,  it  is  difficult  to  determine 
exactly  how  the  belt  forces  interact  with  the  torso  ( i.e., 
where  the  belts  grad)  or  push  on  the  torso).  Figure  3.5 
shows  the  actual  belt  configuration  while  Figure  3.6  shows 
the  final  belt  configuration  used  in  the  simulations  pre- 
sented in  this  report  and  described  below. 

a.  Lap  belt.  The  lap  belt  was  attached  to  the 
anterior- superior  part  of  the  hip  segment  and  anchored  to 
the  sled  behind  and  below  the  point  of  attachment  to  the 
occupant.  The  force  strain  characteristics  were  doubled 
from  that  determined  above  since  there  are  effectively 
two  lengths  of  webbing  restraining  the  subject  (i.e., 
t.he  belt  wraps  around  the  subject  and  is  anchored  at  two 
points  on  the  sled) . 

b.  Upper  and  Middle  Torso  belts.  From  repeated 
observations  of  high  speed  films,  the  impression  was  gained 
that  the  primary  restraining  of  the  torso  is  a result  of 
belt  forces  applied  to  the  chest  in  the  region  of  or  just 
below  the  sternum,  and  that  the  forces  restraining  at  T^^ 
and  the  shoulders  are  not  sufficient  to  prevent  some 
flexion  of  the  torso  above  this  point.  Initially  it  was 
thought  that  these  forces  at  the  chest  were  due  to  the 
chest  belt  which  wraps  around  from  behind  and  passes  just 
beneath  the  arms.  Later  in  the  study.  Dr.  Thomas  indicated 
that  this  belt  is  only  a backup  safety  belt  and  has  too 
much  initial  slack  to  be  of  major  importance.  This  point 
is  academic  to  the  simulation,  however,  since  it  only 
means  that  the  shoulder  belts  have  their  primary  restraining 
action  on  the  chest  rather  than  the  shoulders.  With  these 
points  in  mind,  the  torso  restraint  system  was  modeled  as 
shown  in  Figure  3.6. 
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Figure  3. 


5 Front  and  Side  Photographs  of  NAMRL  Subject  in 
Sled  Chair  Showing  Restraint  System. 
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The  shoulder  belts  are  simulated  by  one  belt  (belt  2) 
which  anchors  slightly  above  and  behind  the  shoulders  and 
attaches  to  the  superior-anterior  region  o£  the  upper 
torso  segment,  and  a second  belt  (belt  3)  which  attaches 
just  below  this  belt  on  the  upper  torso  and  anchors 
vertically  below  near  the  hip.  The  chest  restraining 
forces  are  simulated  by  a rigid  contact  surface^  (line  4 ) 
which  is  specified  such  that  it  is  initially  in  contact 
with  a contact  ellipse  on  the  superior -anterior  part  of 
the  middle  torso  segment,  just  below  joint  3. 

In  order  to  simulate  slipping  of  the  shoulders  rela- 
tive to  the  shoulder  belts,  two  techniques  were  used. 

First,  a model  option  was  implemented  for  belts  2 and  3 
which  limits  the  tension  in  belt  three  to  50%  of  the 
tension  in  belt  2.  Second,  in  order  to  simulate  slipping, 
the  belt  material  properties  of  belt  two  were  altered  from 
the  measured  webbing  properties  for  the  first  5 centimeters 
of  Tj^  movement  as  shown  in  Figure  3.7.  As  illustrated. 


Figure  3.7  Force-Deflection  Specifications  for  Upper 
Torso  Belts. 

; __ — 

Since  the  belt  system  in  the  MVMA-2D  model  is  not 

sufficiently  general  to  represent  the  complex  restraint  sys- 
tem used  in  the  NAMRL  testing,  a fake-belt  was  needed  to 
simulate  the  chest  restraint. 
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the  slipping  force  is  set  equal  to  70  N for  the  belts 
on  the  two  shoulders.  At  5 centimeters  of  slip,  the 
measured  belt  properties  (considering  lengths  and  ntimbers 
of  webbings)  take  effect. 

As  discussed  in  Section  B.l,  the  option  to  include  tor- 
so flexion  resulting  from  belt  slip  at  the  shoulders  is  pro- 
vided for  by  placing  joint  3 just  above  the  chest  restrain- 
ing belt.  For  most  runs  in  this  report,  however,  the  stiff- 
ness coefficient  of  joint  3 is  set  to  a value  which  meUces 
the  torso  essentially  rigid  from  joint  2 to  joint  5. 

c.  Chest  Compliance.  In  comparison  to  the  chest 
compliance,  the  belt  material  may  be  considered  essentially 
rigid  and  so,  for  simplicity,  the  contacting  surface  used 
to  simulate  the  belt  acting  on  the  chest  was  made  rigid. 
Values  of  chest  compliance  in  the  literature  cover  a wide 
range  depending  on  the  condition  of  the  cadaver,  the  rate 
of  force  application,  and  the  size  cuid  mass  of  the  deform- 
ing disk.  Most  experiments  use  a disk  of  ^U^out  6**  dia- 
meter applied  to  the  sternal  area  of  the  chest,  and  it  is 
doubtful  that  the  compliance  factor  arrived  at  in  this 
manner  would  agree  with  that  obtained  by  using  a belt  of 
low  mass  which  contacts  a substantial  surface  of  the  chest. 
As  a result  of  these  considerations,  a value  for  the  chest 
compliance  of  1000  lb. /in.  or  1751  N/cm  v;as  arbitrarily 
specified  for  the  middle  torso  contacting  ellipse.  Effects 
of  this  varying  compli^mce  factor  will  be  illustrated. 

C.  NAMRL  Simulations 

1.  General . This  section  contains  the  graphical 
comparisons  of  the  MVMA-2D  simulations  and  experimental 
results  for  6 and  15  G sled  runs.  Unless  otherwise  noted, 
the  simulations  use  the  data  set  described  in  Section  B of 
this  chapter  developed  from  physical  measurements  on  the 
5 NAMRL  subjects  (described  in  Chapter  2)  and  other  avail- 
able data  and  assumptions.  For  each  run  the  time  dependent 
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variables  of  head  angular  acceleration,  head  angular 
velocity,  head  angular  position,  head  resultant  acceler- 
ation, and  where  appropriate,  Tj  resultant  acceleration 
are  plotted  and  compared  with  experimental  results.  In 
all  cases,  a dashed  line  is  used  for  simulations  while  a 
solid  line  is  used  for  the  average  experimental  curves. 
These  average  curves  are  shown  by  the  dashed  lines  in  each 
plot  of  Appendix  C which  illustrate  the  individual  response 
curves  from  which  these  averages  were  obtained.  Except  for 
the  Ti  resultant  accelerations  at  15  G‘s,  it  is  seen  that 
all  5 siibjects  show  very  consistent  and  similar  responses. 


2 . Results  using  Ti  acceleration  input. 

a.  General.  Initial  attempts  to  simulate  the  NAMRL 
sled  tests  at  6 and  15  G's  were  less  than  satisfactory. 
Because  of  the  uncertainty  in  modeling  the  restraint 
system,  it  was  difficult  to  know  whether  adjustments 
were  needed  in  occupant  parameters  or  in  the  restraint 
system.  In  order  to  separate  these  two  factors  and  allow 
•'tuning''  of  head  and  neck  parameters  for  optimal  matching 
to  sled  test  results,  it  was  decided  to  fix  Tj  rigidly  to 
the  sled  and  use  the  experimentally  determined  linear  Ti 
accelerations  in  the  x and  z directions  as  sled  acceleration 
input  curves.  Use  of  T)  angular  acceleration  was  also 
considered  but  a review  of  the  NAMRL  results  for  this 
signal  and  for  T]  angular  position  suggested  that  this 
was  not  always  a reliable  measurement.  Figures  3.8  and 
3.9  illustrate  these  Tj  acceleration  curves  obtained  by 
averaging  the  data  for  the  five  NAMRL  subjects. 

Figures  3.10  and  3.11  illustrate  the  simulation 
results  for  head  angular  acceleration,  velocity,  and 
position,  and  head  resultant  linear  acceleration  in  compa- 
rison to  the  averaged  sled  test  results  at  6 and  15  G's 
respectively.  These  results  were  obtained  after  some 
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Figure  3.8  Averaged  Experimental  Tj  Accelerations  from  Five  NAMRL  Subjects  at  6 G's 
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Figure  3.10  Siinulation  Results  Using  6 G Tj  Accelerations 
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Figure  3.11  Simulation  Results  Using  15  G T;  Accelerations 
Muscle  Tension  33%  Maximum. 


adjustnents  in  par^uneters  such  as  joint  viscous  friction, 
joint  stop  stiffness,  and  muscle  tension  as  discussed  in 
Section  B of  this  Chapter.  The  initial  neck  angle  is 
70  degrees  from  tiie  horizontal  in  both  cases  and  the  ini-* 
tial  head  angles  are  95.4  and  93.6  degrees  from  the  hori- 
zontal (head  slightly  back)  for  the  6 and  15  G nins  res- 
pectively. From  these  cxirvos,  it  is  seen  that  the  simu- 
lations match  the  experimental  curves  quite  well.  In  both 
simulations,  however,  the  initial  spike  in  head  angular 
acceleration  is  smaller  and  more  rounded  than  in  the  ex- 
perimental results  and  the  head  angular  position  curves 
rise  earlier  but  at  similar  slopes.  At  6 G's  the  angular 
velocity  shows  a plateau  on  the  downward  slope  but  this 
occurs  later  than  in  the  experimental  curve.  While  these 
results  could  perhaps  be  improved  by  further  adjustments 
in  the  pareoneters,  this  was  not  considered  justified  due 
to  t’.e  uncertainty  of  the  signals  themselves  (i.e., 
review  of  the  high-speed  films  indicate  that  the 
accelerometers  have  considerable  movement  relative  to 
especially  in  15  G tests) . 

b.  Effects  of  Huscle  Tension.  The  question  of 
the  importance  of  muscle  mechanics  on  the  dynamic  response 
of  the  head  and  neck  has  been  of  interest  in  recent  years. 
While  this  study  has  not  completely  resolved  the  question, 
some  preliminary  insight  has  been  gained  and  also  the 
validity  of  the  assumption  of  33%  maximum  muscle  tension 
has  been  tested  by  running  the  model  while  changing  only 
the  level  of  muscle  tension.  Figures  3.12  through  3.15 
illustrate  the  results  obtained  if  muscle  tension  is  set 
and  maintained  at  0%  and  100%  of  maximum.  The  general 
characteristics  of  the  curves  are  not  changed  appreciably, 
except  for  head  angular  position. 

At  0%  of  muscle  tension  (i.e.,  no  muscle  input)  the 
angular  positions  are  uncontrolled  and  become  unrealistic, 
while  at  100%  muscle  tension  the  angular  position  is 
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Figure  3.13  Siniulation  Results  Using  6 G Ti  Acceleration 
Muscle  Tension  = lOOi  Maximum. 
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greatly  reduced.  It  is  interesting  that  at  15  G's  the 
effects  of  changing  the  muscle  tension  are  considerably 
less  than  at  6 G's.  AtlSG's/  the  angular  acceleration, 
angular  velocity,  and  head  resultant  acceleration  curves 
are  changed  little  at  100%  muscle  tension  and  fit  the 
experimental  data  extremely  well  at  0%  muscle  tension. 

At  6 G's  the  effects  of  muscle  tension  are  more  dramatic 
in  head  angular  acceleration  and  particularly  head  angular 
velocity. 

These  results  suggest  that  for  the  6 G runs,  the 
muscles  play  a major  role  and  that  the  assumption  of 
33%  muscle  tension  is  reasonable.  At  15  G's,  the  muscular 
effects  are  less  important  and  they  may  in  fact  be  elimi- 
nated by  some  protective  reflex.  The  fact  that  the  angu- 
lar head  position  increases  greatly  at  0%  muscle  tension 
should  not  be  a major  consideration  at  this  time  since  the 
joint  stop  resistance  coefficients  used  probably  do  not 
represent  the  true  physiological  situation.  Further,  the 
condition  of  0%  muscle  tension  is  also  unrealistic  and 
has  been  used  only  to  dramatize  the  effects. 

3 . Results  with  Sled  Acceleration  Input  and  Restraint 
Systenu 

a.  General.  Since  satisfactory  results  were 
obtained  from  the  model  simulations  with  fixed,  it  was 
assumed  that  established  parameter  values  for  the  occupant 
head  and  neck  were  at  least  "in  the  ball  park".  Therefore, 
attempts  at  simulati'^o  the  complete  occupant  with  appro- 
priate torso  and  restraint  system  modeling  were  undertaken. 
The  final  configuration  of  the  restraint  system  used  in  the 
model  is  discussed  in  section  B.7  of  this  chapter  while 
Figure  3.16  shows  the  sled  acceleration  profiles  used 
for  6 and  15  G runs.^ 

^These  profiles  are  for  NAMRL's  high  rate  of  onset, 
long  duration  acceleration  pulses. 
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Figures  3.17  and  3.18  show  the  simulation  results 
and  averaged  experimental  results  at  15  and  6 G's  res- 
pectively using  a constant  33%  mcocimum  muscle  tension, 
a nearly  rigid  torso  (i.e.,  no  bending  at  joint  3)  and 
a chest  compliance  factor  of  1750  N/cm  (1000  Ib/in) . In 
addition  to  the  four  variabl^^ij  shown  previously,  resul- 
tant accelerations  are  also  compeared. 

At  15  G's,  there  is  good  agreement  between  experi- 
mental and  simulated  results.  For  angular  acceleration, 
the  initial  positive  spikes  match  extremely  well  in 
magnitxide  although  in  the  simulation  this  spike  occurs 
about  10  msec  too  soon.  The  initial  negative  spike  is 
of  considerably  larger  magnitude  in  the  simulation.  For 
angular  velocity,  the  magnitudes  of  the  peak  velocities 
match  well  although  the  simulation  curve  reaches  a peak 
about  5 msec  earlier  than  the  experimental  curve  and 
decreases  to  zero  with  a greater  slope,  reaching  zero 
about  25  msec  earlier.  The  angular  position  curves  peak 
at  the  same  time  with  the  simulation  peak  being  about 
.15  radians  or  8.5  degrees  greater.  The  two  curves  rise 
to  their  peaks  with  approximately  the  same  slopes  although 
the  simulation  curve  precedes  the  experimental  curve  by 
about  10  msec.  The  T^^  resultant  curves  match  extremely 
well  considering  the  fact  that  there  is  probably  some 
error  in  the  experimental  curve  (due  to  the  inability  to 
attach  the  accelerometers  rigidly  to  T^^)  and  the  fact  that 
these  curves  are  a result  of  the  complex  interactions  of 
the  restraint  system  with  the  torso  and  dynamic  feedback 
from  the  head.  Both  curves  show  the  bimodal  nature  with 
peaks  of  similar  magnitude  at  similar  times.  The  head 
resultant  accelerations  do  not  match  quite  as  well  but  both 
curves  are  of  a similar  bimodal  nature  matching  well  in 
magnitudes  but  not  as  well  in  times. 
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Figure  3,17  Simulation  Results  for  15  G Sled  Acceleration  - Muscle  Tension 
33^-^  Maximum,  Chest  Compliance  = 1750  N/cm. 
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Figure  3.1B  Simulation  Results  for  6 6 Sled  Acceleration  - Muscle  Tension 
33-  Maximum,  Chest  Compliance  = 1750  N/cm. 
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Figure  3.18  (continued) 


At  6 G's,  the  matching  is  a reasonable  fit,  though  not 
as  close  as  the  IS  G results.  The  simulation  curve  for  an" 
gular  acceleration  shows  a positive  peak  slightly  larger 
and  about  20  msec  earlier  than  the  experimental  curve,  but 
there  is  excellent  agreement  on  the  negative  acceleration 
portion  of  the  curve.  The  angular  velocity  curves  match 
less  well  with  the  simulation  having  a larger  peak  and,  as 
with  the  15  G test,  decreasing  to  zero  at  a greater  rate. 

The  simulation  curve  also  does  not  show  the  plateau  at  about 
150  msec  although  there  is  a slight  change  in  the  slope  at 
this  point.  The  angular  position  curves  match  well  in 
magnitudes,  but  again  the  simulation  curve  peaks  about 
40  msec  sooner  than  the  experimental  curve.  Head  resul- 
tant and  T^  resultant  acceleration  simulation  curves  show 
the  bimodal  characteristics  and  are  of  similar  magnitude 
to  the  experimental  curves  but  the  second  spikes  occur 
earlier  (i.e.,  the  frequency  is  higher)  for  the  simulation 
curves . 

On  the  angular  position  curves,  the  symbols  indicate 
the  tines  at  which  the  subject  has  reached  the  ^oint  stops 
as  determined  by  the  procedures  xn  Section  B.2.  In  both 
cases  the  subject  initially  contacts  the  joint  stop  at  the 
condyle**  rfith  the  head/neck  joint  in  extension  (i.c.  an<ilv  t 
<11  t i ; .ii « • i . J i 1. 1 ? j .» i I ■ . ru’  f ' .js<  >.  < . Af  i . i ’ j;  t s ut  ns* t 
reaches  the  condyle  Toir.r  stop  in  fli'xion  ,it  about  123  msec 
and  then  reaches  the  lower  nock  joint  stoj  in  flexion  at 
about  132  msec.  At  6 G's  the  subject  reaches  the  condyle 
joint  stop  in  flexion  at  about  156  msec  and  does  not  reach 
the  lo«#er  neck  joint  stop. 

b.  Neck  Forces  and  Torques.  Fi;ures  3.19  and 
3.20  s.  jw  the  magnitudes  of  the  necx  'cir.t  tr.rjues  contri- 
1 ited  vy  VISCOUS,  ruscuiar,  and  joint  stop  rec.danisms  ir. 
tne  simulations  at  6 and  15  G's  respectively.  It  wv . 1 bt 
noted  that  the  contributions  due  tc  joint  stops  rv : a*  . 
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Figurt  3.20  N#ck  Joint  Torqut*  in  15  G Simulation.  Muscle 
Tension  • 33i  Maximum*  Viscous  Coefficients  ■ .01  N-m-sec/ 
^ for  Upper  Neck  and  .03  N-m-sec/deg  for  Lower  Neck. 

^int  Stop  Coefficients  - .02^^  N-m/deg  for  Upper  Neck  In 
•<-"»/deg  for  Upper  Neck  In  Extension*  and 
.00S7  N-m/deg  fbr  Lower  Neck  In  Flexion  and  Extension. 
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ly  small  (except  for  the  condyles  in  extension)  du  : to  the 
small  values  of  quadratic  defl 'Action  coefficients  ' sed.  It 
is  also  seen  that  the  muscle  provides  the  major  restraining 
torque  at  both  joints  and  at  both  G levels,  the  proportion 
to  viscous  torques  being  greater  for  the  6 G simulation. 

Figure  3.21  illustrates  the  neck  forces  in  tension 
(negative)  and  the  compression  (positive)  contributed  by 
muscle,  viscous,  and  elastic  components  at  6 and  15  G's. 

It  is  seen  that  the  proportional  contributicjn  of  tiie  force 
due  to  muscle  is  greater  at  6 G's  than  at  15  G's. 

c.  Belt  Forces.  Figur’i  3.22  illuistrates  the  belt 
forces  developed  versus  time  for  each  run.  Since  -he  de- 
flection of  the  upper  torso  belt  does  not  ciitoeed  5 cm 

(it  reaches  about  3.5  cm  at  15  G's)  the  forces  in  >:his 
belt  (belt  2)  and  belt  3 (see  Figure  3.6)  do  not  ex.ceed 
70  N.  As  one  might  expect,  the  force  tiJMe -curve  fur  the 
fake  chest  belt  shows  the  same  multi-peak  characteristics 
as  the  resultant  acceleration  curves. 

d.  Effect  of  Increasing  Neck  Joint  stop  Stiffness. 
Figure  3.23  illustrates  the  effect  on  the  15  G results  of 
xncreasing  the  joint  stop  quadratic  deflection  coefficients 

to  what  would  seem  to  be  more  realistic  values  of  • 

2 ^ 

2.0  N -m/deg  at  the  condyles  and  10  N-m/d  ’ at  C^-T^^. 

The  primary  effect  (compare  with  Figure  3.18)  is  the  large 

negative  spikes  on  the  head  angular  acceleration  curve 

which  result  from  the  head  suddenly  contacting  the  more 

rigid  stops.  It  will  be  noticed,  however,  that  the  peak 

head  angle  position  is  not  reduced  signific2untly  from  that 

which  resulted  from  the  softened  stops  used  in  all  other 

simulations . 

e.  Effect  of  Chest  Compliance.  Figures  3.24 
and  3.25  illustrate  the  effects  of  changing  the  chest 
compliance  factor  to  3500  N/cm  and  875  N/cm  respectively. 

As  expected,  the  primary  change  is  in  th«?  frequency  or 
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Figure  3.?1  (continued) 


3.23  (continued) 


tiJMs  of  occurrence  of  tbe  peeks  and  valleys.  Vith  a 
stiffer  chest  the  curves  are  sharper  and  peaks  occtir 
earlier.  If  one  considers  the  resultant  acceleration 
to  be  the  best  aeasure  of  chest  coopliance,  then  it  would 
aiqpear  that  1751  n/caa  (Figure  3.17}  is  the  best  value.  The 
results  for  875  N/cat  (Figure  3.25)  are  coa^licated»  however, 
by  the  fact  that  this  lower  chest  compliance  allowed 
greater  than  5 cm  deflection  of  the  upper  torso  belt  and 
the  real  belt  properties  cane  into  play  at  about  96  msec. 

(see  Figure  3.7)  This  resulted  in  a sudden  additional 
acceleration  to  the  chest  and  head  causing  the  curves  to 
change  shape  more  than  would  have  been  produced  by  changing 
cheat  compliance  alone.  It  is  clear  that  further  work  is 
needed  to  improve  this  part  of  the  restraint  system  BK>deling. 

f.  Effect  of  Reducing  Condyle  Joint  Stop  Stiff • 
ness  in  Extension.  As  indicated  in  Section  B.3,  the  con- 
dyle joint  stop  stiffness  coefficient  in  extension  was 
maintained  reasonably  close  to  the  MVMA-2D  baseline  values  at 
1.0  N-m/deg^  «diile  other  joint  stop  coefficients  were 
reduced  substantially.  The  reason:  ig  behind  this  was  that 
in  tlie  initial  position  the  subject's  head/neck  angle  is 
very  close  to  the  maximum  head/neck  angle  in  extension 
determined  by  the  procedures  outlined  in  Section  B.2.  Since 
considerable  effort  is  exerted  by  the  subjects  to  achieve 
this  position  during  range-of -root ion  testing,  the  NVKA-2D 
coefficients  seemed  a more  reasonable  approximation  to 
model  the  initial  joint  torque  situation.  In  an  effort 
to  determine  the  appropriateness  of  this  assizmption,  a 
simulation  run  was  made  with  this  joint  stop  coefficient 
reduced  to  .0261  N-m/deg  as  is  used  for  the  condyles  in 
flexion.  Figure  3.26  shows  ttie  results  of  this  run.  In 
comparing  these  curves  with  those  of  Figure  3.18,  it  is 
seen  that  the  only  significant  change  is  ti  at  the  initial 
peak  of  the  angular  acceleration  curve  is  reduced  slightly 
but  more  important,  it  is  delayed  in  time  so  that  it  is 
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Flour®  3.25  SInuUtion  Results  for  15  G Sled  Acceleration- 
»fescle  Tension  ■ 33X  Maximum,  Chest  Compliance  * 875  H/cm. 
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more  in  phase  with  the  experimental  carve.  Figure  3.27 
illustrates  the  head/neck  angle  curves  for  the  t«n>  cases. 
The  difference  is  small  but  significeuit  in  terms  of  the 
time  shift  in  the  peak  of  tl^  angular  acceleration  curve. 
For  the  smaller  joint  stop  stiffness  coefficient,  the  head 
is  allowed  to  extend  backward  relative  to  the  neck  an 
additional  2-3  degrees,  resulting  in  a time  delay  of  about 
8 msec  before  the  head  begins  to  rotate  forward.  The 
effect  that  this  change  would  have  on  improving  the  6 G 
simulations  has  not  yet  been  determined,  but  it  is  evident 
that  further  work  is  needed  in  modeling  the  joint  stop 
characteristics . 
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Figure  3.27  Head/Neck  Angle  versus  Time  for  Condyle  Joint  Stop 
Stiffnesses  of  .0261  N-m/degt  (dashed  line)  and  1.0  N-m/deg? 
(solid  line). 
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g.  Effect  of  Adding  Upper  Torso  Flexion.  As 
described  in  section  B.l,  the  torso  segment  lengths  %fere 
changed  from  the  proportions  of  the  MVNA-2D  baseline  data 
set  so  that  the  upper  torso  joinf  (joint  3 of  Figure  3.1) 
is  above  the  chest  restraint  belt.  In  this  way,  some 
amount  of  torso  flexion  which  is  observed  in  the  high 
speed  movies  may  be  added  by  adjusting  the  stiffness  of 
this  joint.  In  the  runs  presented  so  far,  this  stiffness 
was  maintained  sufficiently  high  (500  N-m/deg)  so  that 
almost  rto  torso  flexion  occurred.  This  was  done  primarily 
because  the  maximum  head  angle  of  the  simulations  was  al- 
ready greater  than  the  experi'  '•ntal  results  and  adding 
torso  flexion  would  increase  this  angle  further.  In 
addition,  the  pareuneters  established  by  using  accelera- 
tions as  the  input  were  based  upon  no  rotation  of  T^.  It 
is  considered  important,  however,  that  this  feature  be 
eventually  included  in  the  model.  Figure  3.28  shows  the 
results  obtained  at  15  G*s  by  reducing  the  linear  stiffness 
coefficient  of  joint  3 from  500  to  75  N-m/deg  while 
maintaining  the  reduced  joint  stop  coefficient  of  the  con- 
dyles in  extension  (Figure  3.26).  Again,  the  changes  are 
not  dramatic  but  are  informative.  The  angular  acceleration 
curve  is  seen  to  be  further  improved  in  that  the  initial 
peak  is  increased  in  magnitude  slightly  and  is  a better 
fit  to  the  experimental  curve  in  the  time  following  this 
peak.  The  angular  velocity  curve  no  longer  has  the 
plateau  at  150  msec  and  crosses  zero  slightly  later  in 
time.  While  the  angular  position  curve  reaches  a greater 
maximum  angle  by  about  6 degrees  (the  upper  torso  angle 
now  flexes  by  about  6 degrees)  it  does  not  sliow  the 
tendency  to  curve  up  again  at  275  msec.  While  these 
changes  represent  improvement  in  the  simulations,  it  is 
curious  that  the  T^  and  head  resultant  acceleration  cux'ves 
match  less  well  to  the  experimental  results  when  this 
torso  flexion  is  allowed. 


Figure  3,28  (continued 
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Chapter  4 

SIMULATIONS  FOR  18>24  YEAR  FEMALES 
A.  18-24  Year  Female  Data  Set. 

While  there  are  several  aspects  of  the  NAMRL  data 
set  and  simlation  %diich  can  be  improved  by  further  work 
and  investigations#  an  attempt  was  made  at  this  time  to 
predict  the  sled  test  results  that  might  be  expected  if 
young  females  %#ere  tested  at  6 and  IS  G's.  The  data 
set  for  the;:e  subjects  was  obtained  by  scaling  the  KANRL 
data  set  parameter  values  using  similar  procedures  to 
those  discussed  in  Chapter  3#  Section  B and  the  measure- 
ment data  obtained  for  these  subjects  and  presented  in 
references  22  and  23. 

1,  Segment  Sgaecificationa 

a.  Torso  and  Extr end  ties.  Torso  length  was 
computed  in  the  same  manner  as  for  the  NAMRL  subjects 
and  the  segment  lengths  taken  in  the  same  percentages 
of  total  torso  length.  lengths  were  taken 

from  the  same  anthropomc'tx  -c  rwsiisurements  and  again 
the  arms  were  removed  at.d  a i-s  of  the  upper  arm 

mass  was  added  to  the  upper  mas^r.  Values  of 

mass#  moment  of  inertia,,  and  dist4.\‘^u  to  centers  of 
mass  from  joints  were  computed  by  sci^Iiiig  to  the 
NAMRL  data  set  values  in  a manner  simlli..r  to  that 
used  for  scaling  NAMRL  data  to  the  baseline  data. 

The  result  of  the  calculations  gave  the  parameter 
values  shown  in  Table  4.1. 


Table  4.1 


Torso  and  Extremity  Segment  Specifications 
For  18-24  Year  Female  Data  Set 

End  of  Link  to  2 

Segment  Length (<r^)  Center  of  Mass  (cm)  Masa(kq)  lYY (Itg-ei'^) 


Upper  Torso 

17.78 

8.89 

11.4 

.087 

Middle  Torso 

22.95 

8.5 

7.9 

.095 

Hip 

9.95 

3.84 

6.55 

.133 

Upper  Leg 

41.88 

19.03 

13.7 

.222 

Lower  Leg 

40.5 

25.47 

7.4 

.208 

b.  Head  and  Neck  Mass  and  Monent  of  Inertia. 

Head  mass  for  the  females  was  estimated  by  assuming  that 
head  mass  is  proportional  to  [Head  Length  x Head  Breadth] 
and  scaling  to  the  Navy  data  based  on  the  ratio  of  there 
valxies  for  both  data  sets.  By  these  calculations,  the 
female  head  mass  is  .885  times  the  NAMRL  head  mass  giving 
an  instrumented  head  mass  of  4.164  kg  for  the  females. 

As  with  the  NAMRL  subjects,  female  head  moment  of  inertia 

was  calculated  from  the  male  cadaver  data  of  Chandler,  et 

al.  (3)  using  the  relationship  that  head  moment  of  inertia 

2 2 

is  proportional  to  [(Menton  to  vertex)  (head  length)  1 

X [Head  circumference] . This  gives  an  instrumented  head 

2 

moment  of  inertia  for  young  females  of  .0211  Kg-m  . 

Female  neck  mass  was  scaled  to  the  NAMRL  neck  mass 
by  assuming  that  mass  is  proportional  to  volume  and  that 
the  neck  is  a cylinder  whose  circumference  is  equal  to 
the  average  of  superior  and  inferior  ncKSk  circumference  and 
%d\ose  length  is  proportional  to  erect  sitting  height. 

This  gives  a female  neck  mass  of  1.194  Kg.  of  which  334 
was  placed  at  the  condyles  and  67%  at  C^-T^. 

c.  NcK:k  Length  and  Location  of  Head  Center  of 
Gravity.  The  location  of  the  head  center  of  gravity  rela- 
tive to  the  condyles  was  determined  by  measuring  distances 


106 


from  tragion  to  the  condyles  on  x-rays  from  young  females 
in  the  IIHS  sagittal  study.  This  gave  average  distances 
for  tragion  of  .987  cm  forwaurd  aj)d  1.966  cm  above  the  con- 
dyles. The  distances  given  by  Ewing r et  al.  (8)  for  the 
location  of  the  c.g.  relative  to  tragion  and  the  shift  in 
the  c.g.  due  to  the  instruaient  package  (9)  on  young  males 
were  assumed  for  the  females  also,  and  were  add<K3  to  the 
distances  of  tragion  to  condyles  to  give  the  location  of 
the  head  center  of  gravity  relative  to  the  coriyles.  This 
gave?- 

X distance  condyles  to  head  c.g.  » 2.3  cm. 

y distance  condyles  to  head  c.g.  ■ -4.07  cm. 

Neck  lengths  for  young  females  were  obtained  by 
direct  measurement  on  x-rays  of  the  linear  distance  from 
the  condyle  to  ^7"^!  disk.  After  appropriate  scaling 
this  give  an  average  neck  length  of  10.6  cm.  This  is  in 
good  agreement  with  the  value  of  10.4  cm  which  would  be  ob- 
tained if  neck  length  were  scaled  to  the  NAMRL  neck  length 
by  the  proportions  of  erect  sitting  height  [ i . e . , 

(female  sit.  ht. /NAMRL  sit.  ht.)  x NAMRL  neck  length  - 
(86.24/92.88)  x 11.2  - 10.4  cm). 

2.  Head  and  Neck  Range  of  Motion.  Values  of  a,  g, 

Y , f',  were  determined  for  the  NAMRL  data  set  from  x-rays 

of  young  males  in  the  IIHS  sagittal  study  (see  Section  B.2). 
The  mean  ranges  of  motion  in  flexion  and  extension  for 
this  group  were  62.5  and  79.6  degrees  respectively.  For 
the  young  females  the  mean  values  for  flexion  and  extension 
were  60.9  and  77.1  degrees  respectively.  Since  these 
values  do  not  differ  significantly  from  the  values  for  the 
males,  the  s^tme  values  for  a,  g,  y,  axid  6 were  used  for  the 
female  data  set. 

3.  Passive  Joint  Torques  and  Joint  Stops.  Since  the 
joint  stop  quadratic  deflection  coefficients  used  for  the 


NAMRl.  su!i>octs  havo  physiolcxjica  1 basis,  the  soMC 

values  uare  used  for  the  young  fenales.  Viscous  co- 
effiricats,  ho%#ever,  were  scaled  to  the  NAMRL  data  set 
valuwss  of  .01  N-m-sec/deg  and  .03  N-m-sec/deg  for  the 
upper  and  loifer  neck  joints  respectively  by  the  propor- 
tion of  neck  cross-sectional  areas.  The  average  of  sixpmriox 
and  inferior  neck  circtwferences  vas  used  to  detaureine  neck 
radius  which  was  cooputed  to  be  5.99  cm  for  the  NAMBL 
subjects  and  5.49  cm  for  the  fenales.  This  gave  a scaling 
factor  of  .84  which  gives  values  of  .0084  and  .0252  N-n- 
sec/deg  for  the  female  upper  and  lower  neck  viscous  co- 
efficients. 

4.  Neck  Muscle  Parameters.  Neck  muscle  parameters 

i I j, 

ai,  a^,  S3,  and  ai,  a;,  aa  were  scaled  to  the  NAMRL  para- 
meters using  the  relations  given  in  equations  1 through  4 
of  section  B.4  of  Chapter  3.  Distances  l2$  and 

(see  Figure  3.4)  %icre  assumed  to  be  the  same  for  the 
young  females  as  were  estimated  for  the  NAMRL  subjects. 

The  neck  length  of  10.6  cm  was  used  for  In  extension, 

the  average  pull  force  for  young  females  was  27.04  Ibf 
compared  to  48.3  Ibf  for  the  NAMRL  subjects.  In  flexion 
the  average  pull  force  for  young  females  was  19.4  Ibf 
compared  to  36  Ibf  for  NAMRL  subjects.  Using  these 
values  and  the  procedures  outlined  in  section  B,4.  of 
Chapter  3 and  in  greater  detail  in  reference  1,  the  muscle 
parameter  values  shoim  in  Table  4.2  for  the  young  fesule 
data  set  were  calculated. 
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Table  4.2 
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Muscle  Paraneter  Values  Used  in  18-24  Year  Female  Data  Set 


Par2uneter 


100%  Muscle 

Location  aj^(N-m/deg)  ajCdeg”^)  a^  (sec/deg)  Torque  (N-m) 


Occ.  Condyles  .03 

.086 

.007 

6.51 

C^-Tj^  .093 

.086 

.007 

19.99 

a^  (N/cm) 

0 

a^ (sec/cm) 

100%  Muscle 
Tension  (N) 

Neck  Stretch  73.7 

1.29 

.15 

938.5 

H 


0^ 


% 

'IfT' 


t 


As  with  the  NAMRL  subiect.s  the  muscle  torques  and 
tensions  <were  maintained  constant  at  33%  of  maximum  through- 
out the  simulati'>n . 

5 . Neck  Stretch  and  Compression  Parameters.  The 
elastic  and  viscous  coefficients  which  describe  the  visco- 
elastic properties  of  the  neck  in  stretch  and  compression 

were  maintained  at  C_  » 10.55  N-ssc/cm  and  «■  1636  N/cm 

s s 

for  the  young  female  subjects.  Variations  of  these  values 
in  simulations  of  NAHRL  subjects  showed  the  itiodel  to  be 
relatively  insensitive  to  changes  in  these  pareuneters. 

6.  Head  and  N€Kik  Initial  Angles.  The  initial 
position  of  the  head  and  neck  were  maintained  the  same  as 
for  the  NAMRL  subjects. 

7.  Restraint  Syst<g«  and  Chest  Compliance.  The 
restraint  system  configuration  was  maintained  the  same  as 
for  NAMRL  subjects  except  for  reposi (:ioning  the  anchor 
points  and  attachment  points  due  to  the  change  in  body 
segment  sizes.  Chest  compliemcc  was  maintedned  at  1751  M/cm 
or  1000  Ib/in. 
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B.  Simulations  for  18-24  Year  Females. 


Figures  4.1  and  4.2  show  the  simulation  results 
(dashed  lines)  in  comparison  to  the  NAHRL  experimental 
results  :olid  line)  and  NAMRL  simulation  results  (dotted 
line)  for  6 and  15  G sled  tests.  In  general,  the  results 
are  not  markedly  different  from  the  NAMRL  results.  The 
primary  difference  is  in  the  angular  position  curves  where 
the  maximum  head  rotation  is  approximately  40  and  25  per- 
cent greater  in  the  6 and  15  G runs  respectively.  This 
increase  in  peak  head  angle  is  probably  largely  due  to 
the  weaker  neck  muscles  and  consequent  changes  in  neck 
muscle  model  pareuneters.  A second  probable  consequence 
of  the  weaker  muscles  is  seen  in  the  angular  velocity 
curve  at  6 G's  where  the  curve  is  seen  to  return  to  zero 
more  slowly  (and  interestingly  in  closer  approximation 
with  the  NAMRL  experimental  cuirve) » At  15  G's  this  portion 
cf  the  angular  velocity  curve  is  not  changed  as  much  and 
this  is  perhaps  another  indication  that  the  role  of  neck 
muscles  is  less  important  at  higher  G levels. 

At  6 G's  the  initial  peak  in  the  angular  acceleration 
curve  is  nearly  identical  in  magnitude  and  time  to  tne 
NAMRL  simulation,  but  both  occur  sooner  in  time  than  the 
experimental  results.  At  15  G's  the  initial  peak  for  the 
females  is  sligh  .ly  larger  and  occurs  slightly  earlier 
than  the  NAMRL  simulation.  This  earlier  peaking  may  be  a 
consequence  of  the  smaller  head  mass  and  moment  of  inertia 
for  the  females.  At  both  6 and  15  G's  the  initial  negative 
peak  is  reduced  from  the  NAMRL  simulations  and  this  may 
also  be  a consequence  of  the  reduced  muscle  strength.  The 
reduction  in  these  negative  peaks  is  seen  to  be  greater 
at  6 G's  than  at  15  G's. 

Concerning  the  resultant  acceleration,  it  is  seen  that 
at  6 G's  the  frequency  of  the  bimodal  portion  of  the  curves 
is  increased  from  the  NAMRL  simulations  although  the 
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initial  spike  occurs  at  about  the  same  time.  At  15  G's 
the  initial  spike  occurs  sooner  but  the  frequency  is 
not  significantly  different. 


Ill 


Figure  4.1  (continued) 
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Figure  4.1  (continued) 
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Figure  4.2  15  G Simulation  Results  for  18-24  Year  Females  (dashed  lines)  in 
Comparison  with  NAHRL  Experimental  (solid  lines)  and  NAHRL  Simulation  (dotted  lines). 


Figure  4.2  (continued) 
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Chapter  5 

DISCUSSION  AND  RECOKHENDATIONS 

Ideally,  the  fulfilltsient  cf  the  objectives  ot  tliis 
research  effort  will  ultinately  result  in  1}  a better 
understardiiag  of  the  r lationshipe  and  inportance  of 
various  physical  characteristics  to  the  dynamic  respossse 
(and  therefore  injury  susceptibility)  of  the  head  and  neck, 
and  2)  a mathesuitical  model  which  utilizes  thorn  parameters 
and  which  will  enable  accurate  and  reliable  predictions  of 
head  and  neck  responses  over  a large  rajuge  of  conditions 
aund  variations  in  subject  physical  characteristics.  The 
results  and  accomplishments  of  the  initial  12  months  of 
this  study  have  been  satisfying  and  encouraging  toward 
achieving  these  goals. 

The  .iimulation  results  of  Chapter  3 illustrate 
excellent  reproductions  of  NAMRL  sled  test  results  for 
head  angular  acceleration,  head  angular  velocity,  head 
angular  position,  and  head  and  Ti  resultant  acceleration 
curves.  Indeed,  it  might  be  concluded  that  these  results 
are  sufficient  and  that  the  model  can  now  be  used  to  obtain 
reasonable  predictions  for  tasks  6 through  8 (see  section  C, 
Chapter  1) . It  is  believed,  however,  that  there  are  sig- 
nificant improvements  which  can  still  be  made  in  the  model 
validation  which  will  both  add  to  our  understanding  of  the 
mechanisms  involved  in  a dynamic  situation  and  improve  the 
capability  and  credibility  of  the  model. 

One  area  of  particular  concern  in  the  simulations 
presented  to  date  is  in  the  modeling  of  joint  stops  and 
passive  tissue  resistance,  and  the  manner  in  which  measured 
voluntary  range-of -root ion  limits  relate  to  these  phanonsma. 
While  the  simulations  presented  in  Chapter  3 are  reasonable 
fits  to  experimental  results,  there  is  considerable 
question  as  to  the  values  of  joint  stop  parameters  used. 
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Hiis  question  nust  be  resolved  if  tlie  model  is  to  be  a 
relii&ble  predictor  for  persons  with  smaller  ranges  of 
motion  (e.g.,  elderly  persons)  or  for  higher  G levels 
where  the  stop  plays  a more  significant  role  in  res- 
training the  head.  A more  complete  resolution  of  the 
problem  will  require  experimental  measurements  on  husuuis 
and  animals / but  it  is  also  likely  that  further  attesq>ts 
at  simulating  these  charact<?ri»tics  will  lead  to  more 
suitaUale  and  realistic  simulation  results  and  a better 
understanding  of  this  area. 

The  importance  of  muscle  forces  in  the  dynamic 
response  of  the  human  neck  has  been  investigated  to  some 
degree  in  this  study,  but  there  is  more  to  be  done. 

Results  so  far  suggest  that  muscle  effects  xre  more  sig- 
nificant at  6 G's  than  at  15  G's  on  the  NAMRL  sled  tests 
but  that  subjects  may  not  be  using  their  muscles  to  the 
maximum  extent  possible  even  during  the  15  G sled  runs. 

While  an  increase  in  muscle  tension  above  the  33%  of 
maximum  used  would  improve  the  simulation  head  angular 
position  curve  by  reducing  the  peak,  it  would  also  result 
in  a more  rapid  decline  of  head  angular  velocity  which 
is  already  too  steep.  Further  work  with  this  aspect  of  the 
model  should  prove  useful  in  understanding  the  role  of 
muscle  and  improving  the  model  performance.  The  use  of 
EMG  signals  during  dynamic  testing  and  experiments  with 
animals  would  provide  useful  information  for  understanding 
the  action  of  muscle  under  dynamic  conditions. 

For  the  NAMRL  sled  test  simulations  and  predictions, 
further  work  with  the  restraint  system  is  needed,  especially 
with  regard  to  the  apparent  slipping  of  the  shoulders 
relative  to  the  shoulder  belts  and  the  consequent  torso 
bending.  Measurements  of  belt  forces  during  sled  tests 
would  be  extremely  useful  to  improving  this  aspect  of  the 
model  which  is  essential  if  sled  test  results  are  to  be 
accurately  extended  to  the  general  adult  population. 
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In  Chapter  1,  section  C,  eight  tasks  %fere  listed  which 

would  lead  to  achievenent  of  the  project  goals.  Of  these « 

the  first  three  have  been  coupleted  and  a substantial  aai 

prcsoising  start  ha:,  been  made  on  task  four.  In  view  of  these 

<» 

achievwssnts  and  the  results  to  date*  it  is  recosmended  that 
efforts  toward  cosipleting  task  four  (as  outlined  in  the  above 
discussion)  be  continued*  and  where  appropriate  and  feasible 
that  experiments  be  conducted  to  assist  with  the  OK>del  vali** 
dation.  As  an  additional  part  of  this  validation*  it  is 
also  recomnended  that  efforts  be  included  to: 

1)  extend  simulations  to  include  3*  10*  and  greater 
than  15  G (if  avail2d>le)  sled  tests. 

2)  include  simulations  of  sled  tests  with  accelera- 
tion profiles  of  different  rates  of  onset  and 
different  durations. 

3}  include  simulations  of  sled  tests  with  different 
head/neck  initial  positions. 

At  the  sasuB  time  work  on  task  five  can  be  started.  This 
will  involve : 

1}  studying  experimental  results  of  NANHL  subjects 
vdiose  physical  characteristics  differ  from  the 
group  used  for  simulations  of  this  report. 

2)  extending  simulations  to  these  other  NANRL 
subjects. 

3}  performing  statistical  correlations  of  physical 
measurement  results  with  peak  parameter  values  of 
sled  tests  for  the  18  NANRL  subjects  measurec!. 

Upon  completion  of  these  tasks*  work  on  tasks  six 
through  eight  can  be  xindertaken  with  the  expectation  that 
the  predictions  will  be  reliable  extrapolations.  The  results 
will  therefore  be  useful  toward  defining  the  envelopes  of 
impact  acceleration  which  result  in  injury.  In  addition* 
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by  these  procedures  vdiich  conbine  physical  neasureaents,  ex' 
perimental  data,  and  matheiaatical  modeling#  an  increased 
understanding  of  the  factors  influencing  dynamic  responses 
during  impact  can  be  achieved.  This  will  provide  iaqportant 
information  for  design  of  improved  dummies  for  impact 
investigations . 
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APPENDIX  A 


CROSS  REFERENCE  TABLES  FOR  MEASUREMENT  CODE  NAMES 

Tables  A.l  through  A. 4 provide  a cross  reference  for 
the  abbreviated  measurement  names  used  in  the  t2U3les  of 
statistical  results  in  Chapter  2 and  in  the  tables  of 
measurement  results  by  subject  in  Appendix  B. 


'u' 

[ 

,! 

TABLE  A.l 
CODE  NAME 

ANTHROPOMETRY  CODE  NAME  CROSS  REFERENCE 
MEASUREMENT 

1 

WT(KR) 

Weight  in  kg 

1 

MT(LB) 

Weight  in  lb 

1 

STAT(CM) 

Stature 

1 

PONDINDX 

Ponderal  Index 

1 

ERSITHT 

Erect  Sitting  Height 

) 

HEADCIR 

Head  Circumference 

HEADELPS 

Bennett  Ellipse  Ci rcumference 

BITRfiOI 

Bitraqion  Diameter 

HEADBR 

Head  Breadth 

HEADLG 

Head  Length 

SAGARC 

Sagittal  Arc  Length 

CORARC 

Coronal  Arc  Length 

1 

BITRGGLB 

Bi tragi on-Glabel la  Arc  Length 

BITRGMEN 

Bi tragi on-Nenton  Arc  Length 

6ITRGINA 

Bitragion-Inion  Arc  Length 

FACEHT 

Facial  Height 

'1 

LATNKBR 

Lateral  Neck  Breadth 

V 

APNKBR 

Anterior-Posterior  Neck  Breadth 

V 

SUPNKCIR 

Superior  Neck  Circumference 

' 

IMFNKCIR 

Inferior  Neck  Circumference 

P0STNKL6 

Posterior  Neck  Length 

BIACRBR 

Bi acromial  Breadth 

BIOELTBR 

Shoulder  Breadth  (Bi deltoid} 

■'> 

CHESTHT 

Chest  Height 

CHESTBR 

Chest  Breadth 

CHESTCIR 

Chest  Circumference 

WAISTHT 

Waist  Height 

WAISTBR 

Waist  Breadth 

WAISTCIR 

Waist  Circumference 

HIPHT 

Hip  Height 

HIPBRSTD 

Hip  Breadth  (Standing  Erect) 

HIPCIR 

Hip  Circumference 

ACCRRADLG 

Acromion-Radiale  Length 

ARMCIRAX 

Upper  Arm  Circ.  (at  Axilla) 

ARMCIREL 

Upper  Arm  Circ.  (above  Elbow) 

BICFLCIR 

Biceps  Flexed  Circumference 

i 

RADSTYLG 

Radiale-Stylion  Length 

1 

FRARMCIR 

Forearm  Circumference 

1 

WRISTCIR 

Wrist  Circumference 

( 

,■ 

HANOLG 

Hand  Length 

I 

TRCFEMLG 

Trochanter-Femoral  Condyle  Length 

1 

UPTHKIR 

Upper  Thigh  Circumference 

! 

”■ 

LWTHICIR 

Lower  Thigh  Circumference 

; i 

FIBULALG 

Fibula  Length 

il 

i' 

FIBULAHT 

Fibula  Height 

i 

■“:  ’ 
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TABLE  A.1  (continued) 


CODE  NAME 


MEASURENENT 


CALFCrR 

ANKLECIR 

FOOTLG 

FOOTER 

HUMOIA 

FEMDIA 

TRICPSF 

SUBSCPSF 

SUPILSF 

LTIL-SYM 

RTIL-SYM 

ASISBR 

NRMSITHT 

TRAGHTS 

TRAGDPS 

GLABHTS 

BLA30PS 

EYELPHTS 

EYELPDPS 

EYELPWOG 

C7HTS 

C7DPS 

SSTRNHTS 

SSTRNDPS 

SHLDRHTS 

SHLOROPS 

SHLDRBR 

ILCSPHTS 

ILCSPDPS 

BISPNBR 

TRCHHTS 

TRCHOPS 

BITRCHDI 

HIPBRSIT 

ORBHn 

ORBDPT 

TRAGHTC7 

TRAGOPC7 

GLABHTT 

aABDPT 

EYELPHTT 

EYELPDPT 

ECTCNATT 

ECTCNDPT 


Calf  Ci rcumference 
Ankle  Circumference 
Foot  Length 
Ball-of-Foot  Breadth 
Humeral  Biepi condylar  Oia. 

Femoral  Biepi condylar  Oia. 

Triceps  Skinfold  (mm) 

Subscapular  Skinfold  (mm) 

Suprailiac  Skinfold  (nai) 

Left  Asis  to  Symphyslon 

Right  Asis  to  Syrophysion 

Anterior  Superior  Iliac  Spine  (Asis)  Breadth 

Normal  Sitting  Height  (re  SRP) 

Tragion  Height  (re  SRP) 

Tragion  Depth  (re  SRP) 

Glabella  Height  (re  SRP) 

Glabella  Depth  (re  SRP) 

Eye  Ellipse  Point  Height  (re  SRP) 

Eye  Ellipse  Point  Depth  (re  SRP) 

Eye  Ellipse  Point  Width  (re  Glabella) 
Cervicale  Height  (re  SRP) 

Cervicale  Depth  (re  SRP) 

Suprasternal e Height  (re  SRP) 

Suprasternale  Depth  (re  SRP) 

Shoulder  Height  (rf  SRP) 

Shoulder  Depth  (re  SRP) 

Shoulder  Breadth 

Anterior  Superior  Iliac  Spine  Ht  (re  SRP) 
Anterior  Superior  Iliac  Spine  Depth  (re  SRP) 
Bi spinous  Breadth 
Trochanter  Height  (re  SRP) 

Trochanter  Depth  (re  SRP) 

Bi trochanter  Diameter 
Hip  Breadth  (Seated  Erect) 

Infraorbitrale  Height  (re  Tragion) 
Infraorbitrale  Depth  (re  Tragion) 

Tragion  Height  (re  Cervicale) 

Tragion  Depth  (re  Cervicale) 

Glabella  Height  (re  Tragion) 

Glabella  Depth  (re  Tragion) 

Eye  Ellipse  Point  Ht  (re  Tragion) 

Eye  Ellipse  Point  Depth  (re  Tragion) 
Ectocanthus  Height  (re  Tragion) 

Ectocanthus  Depth  (re  Tragion) 
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CODE  MME MEASUREMENT 


■iT; 

SHLDRSX 

Shoulder  Point 

~X  Direction 

SHLDRSY 

Shoulder  Point 

-Y  Direction 

SHLDRSZ 

Shoulder  Point 

-Z  Direction 

, "F 
.'.f  ' 

C7 

SX 

Cervicale 

-Y  Direction 

' 

i 

C7 

SY 

Cervicale 

-Y  Direction 

C7 

SZ 

Cervicale 

-Z  Direction 

SSTRNSX 

Suprasternal  e 

-X  Direction 

*V  •'Ifl 

SSTRNSY 

Suprasternal e 

-Y  Direction 

SSTRNSZ 

Suprastemale 

-Z  Direction 

TRAG  SX 

Tragi on 

-X  Direction 

TRAG  SY 

Tragi on 

-Y  Direction 

i'  / 

Trag  SZ 

Traglon 

-Z  Direction 

1%?. 

ORBITSX 

Infraorbitale 

-X  Direction 

ORBITSY 

Infraorbitale 

-Y  Direction 

ORBITSZ 

Infraorbitale 

-Z  Direction 

■vj'- 

GUB  SX 

Glabella 

-X  Direction 

GLAB  SY 

Glabella 

-Y  Direction 

V#:: 

>])«« 

GLAB  SZ 

Glabella 

-Z  Direction 

1 ') 

EYELPSX 

Eye  Ellipse  Point 

-X  Direction 

EYELPSY 

Eye  Ellipse  Point 

-Y  Direction 

EYELPSZ 

Eye  Ellipse  Point 

-Z  Direction 

..tf,, ' 

ECCANSX 

Ectocanthus 

-X  Direction 

K.t 

ECCANSY 

Ectocanthus 

-Y  Direction 

ECCANSZ 

Ectocanthus 

-Z  Direction 

It 

rA  ■ 
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TAfiLE  A. 3 RANft£-QF-$#9TI0N  CODE  WANE  CROSS  REFERENCE 


CODE  mu 


HEASUREMENT 


P2i\'EUTY 

P2NEUTO 

P2NEUTR 

P3EXTY 

P3EXTP 

P3EXTR 

P4FLEXY 

P4FLEXP 

P4aEXR 

P5RTR0TY 

P5RTR0TP 

P5RTRQTR 

P6LTR0TY 

P6LTR0TP 

P6LTR0TR 

P7RLBN0Y 

P7RLBN0P 

P7RLBNDR 

P8LLBN0Y 

PaiBNOP 

P8LLBN0R 

P9LR0aV 

P9LR0FLP 

P9LR0FLR 

PIOLROBY 

PIOLROBP 

P10LR0BR 

PIT  PROXY 

PllRROXP 

P11RROXR 

PSAGROM 

PROTROM 

PLATROM 


ir>)\'oto  2->Naitra1  Head  Position— Yam 
Photo  2— Naitrai  Head  Position— Pitch 
Photo  2— Neutral  Head  Position— Roll 
Piioto  3— Extension— Yaw 
Photo  3— Extension— Pitch 
Photo  3— Extension— Roll 
PixJto  4— Flexion— Yaw 
Photo  4— Flexion— Pitch 
Photo  4— Flexion>Ro17 
Photo  5— Right  Rotation— Yaw 
Photo  5— Right  Rotation— Pitch 
Photo  5— Right  Rotation— Roll 
Photo  6— Left  Rotation— Yaw 
Photo  6— Left  Rotation— Pitch 
Photo  6— Left  Rotation— Roll 
Photo  7— Right  Lateral  Bend— Yaw 
Photo  7— Right  Lateral  Bend— Pitch 
Photo  7— Right  Lateral  Bend— Roll 
Photo  8— Left  Lateral  Bend-Yaw 
Photo  8— Left  Lateral  Bend— Pitch 
Photo  8— Left  Lateral  Bend— Roll 
Photo  9— Left  Rotation  + Flexion— Yaw 
Photo  9— Left  Rotation  Flexion— Pitch 
Photo  9— Left  Rotation  * Flexion— Roll 
Photo  10— Left  Rotation  + Left  Lateral  Bend— Yaw 
Photo  10— Left  Rotation  + Left  Lateral  Bend— Pitch 
Photo  10— Left  Rotation  + Left  Lateral  Bend— Roll 
Photo  11— Right  Rotation  + Extension— Yaw 
Photo  11— Right  Rotation  + Extension— Pitch 
Photo  11— Right  Rotation  + Extension— Roll 
Sagittal  Range  of  Notion  from  Photogranmetry 
(P3EXTP  + P4FLEXP) 

Rotational  Range  of  Notion  from  Photogranmetry 
(P5RTR0TY  + P6LTR0TY) 

Lateral  Bend  Range  of  Motion  from  Photogranmetry 
(P7RLBNDR  + P8LLBN0R) 
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TABLE  A.4  REaEX  TIMES  AMD  STRENGTH  CODE  NAME  CROSS  REFERENCE 
CODE  NAME  MEASURD€N7 


RFL  LAT 
RFL  PLXR 
RFL  EXTR 
RFL  AV 
STR  RTL 
STR  LTL 
STRUTAV 
STR  EXTR 
STR  FLXR 
STRSAGAV 
SAGLATAV 


Reflex  Tiae  In  Lateral  Flexion 

Reflex  Tine  of  Flexor  Muscles 

Reflex  Time  of  Extensor  Muscles 

Average  of  All  Reflex  Tines 

Pull  Force  Fran  Right  Lateral  Flexors 

Pull  Force  Fron  Left  Lateral  Flexors 

Average  Pull  Force  Fron  Left  and  Right  Lateral  Flexors 

Pull  Force  Fron  Extensor  Muscles 

Pull  Force  Fron  Flexor  Muscles 

Average  Pull  Force  Fron  Extensors  and  Flexors 

Average  Pull  Force  Fron  Lateral  and  Sagittal  Tests 
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APPENDIX  B 

MEASUREMENT  RESULTS  BY  SUBJECT 

Tabltttf  B.I  through  B.4  give  the  meaeureaent  resulta 
by  individual  subject  for  anthropometry head  ^md  torso 
landmarks/  range  of  motion/  and  strength  and  reflex  time 
results  respectively.  The  «U>breviated  measurement  names 
may  be  cross  referenced  in  the  Tables  of  Appendix  A. 
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TABLE  B.l  (continued) 


TABLE  B.l  (continued) 


TABLE  B.l  (continued) 


TABLE  B.2  THREE  SPACE  LOCATIONS  07  UPPER  TORSO  AMD  HEAD 
LAMDMARKS  BT  SUBJECT 
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TABLS  B.3  RAKGE-QF«HOTIOir  AXGLBS  BT  SUBJECT 
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APPENDIX  C 


EXPERIMENTAL  SLED  TEST  RESULTS 


Figures  C.l  through  C.5  show  the  6 G experimental 
time  traces  for  head  angular  acceleration/  heed  angular 
velocity,  head  angular  position,  head  resultant  accelera- 
tion, and  T^  resultant  acceleration  for  the  five  NAMRL 
subjects  used  in  the  simulations  of  Chapter  3.  The  dashed 
line  in  each  figure  is  the  average  of  the  five  curves  and 
corresponds  to  the  solid  line  shown  in  the  simulation  com- 
parisons of  Chapters  4 and  5.  Figures  C.6  through  C.IO 
show  similar  experimental  curves  for  15  G sled  tests. 
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Figure  C,1  6 G Experiment 

Five  NAMRL  Subjects.  Dash 
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Figure  C.9  15  G Experimental  Head  Resultant  Acceleration  Curves  for 
5 Five  NAMRL  Subjects.  Dashed  Line  is  Average  Curve. 
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